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SUMMARY 


This final report presents the results of a study conducted by the 
A (Research Manufacturing Company of California of the preliminary design of 
an alternate heat source assembly (HSA) intended for use in the Brayton 
isotope power system (BIPS) presently under development by the AiResearch 
Manufacturing Company of Arizona under a Department of Energy contract. 

The BIPS converts thermal energy emitted by a radioactive heat source into 
electrical energy by means of a closed Brayton cycle. The HSA, a major 
component of the BIPS, accomplishes the transfer of thermal energy from the 
heat source to the Brayton cycle working fluid. 

The study effort involved the selection and optimization of a heat source 
heat exchanger (HSHX) configuration and material with enhanced reliability and 
fabr icabi I ity as compared with the original concept. The thermal performance 
of the entire HSA was analyzed for various modes of system operation. Particu- 
lar emphasis was placed on the predicted temperature of the iridium post- impact 
containment shell (PICS), the protective housing for the radioactive fuel. 
Required PICS temperatures were specified for the various system operating 
modes. Structural characteristics, including low-cycle fatigue life, natural 
frequencies, and dynamic loading, were analyzed for the HSA. Recommendations 
are made for future analytical and experimental work. A layout drawing was 
prepared for the HSA. 

The study identified a Hastel loy X single- tube helical coil as the HSHX 
configuration that maximizes reliability and fabr icabi I ity within the frame- 
work of the BIPS requirements. The optimized design, which met performance 
requirements, was a 10-turn coil made of a tube of 3.96-cm (1.56-in.) outside 
diameter and 3.703-cm (1.458-in.) inside diameter. 

For the auxiliary cooling system, required during nonoperational phases 
of the BIPS, neon is the recommended fill gas. Transient startup conditions 
have been analyzed, and a low-cycle fatigue life of 570 cycles (well above 
the required 100 cycles) has been predicted for the HSHX. 

PICS temperature requirements can essentially be met for all modes of 
operation except for the emergency cooling system, which is required to cool 
the heat source after a system failure. Analysis indicates that an insuffi- 
cient number of insulation foils are melted to lower the PICS temperatures 
suf f ic i entl y. 

With some minor modifications to the existing HSA structure, natural 
frequency and dynamic loading requ i rements can be met. 
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PROGRAM DESCRIPTION AND OBJECTIVES 
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TABLE I 

DOE-SPECIFIED PICS OPERATING TEMPERATURES 


Operating mode 

Maximum, *C 

Minimum, °C 

Ground hand 1 i ng 

1200 

None 

Launch and ascent (auxiliary cooling) 

1200 

1 050 

Normal operation 

1200 

mo 

Emergency cooling transient 

1600 

1 1 60 

Emergency cooling steady state 

- — 

1200 

1 1 50 


The steady-state lower PICS temperature limit of n«snv k * 

and this requirement should probably be reviewed. restrictive 


The study program was designed to meet the following objectives: 

(1) Selection of an alternate HSHX configuration and material that would 
provide enhanced reliability and fabr I cab 1 1 ity as compared with the 
former concept 

(2) Optimization of the selected design with respt." to weight, pressure 
drop, and working fluid outlet temperature 

(3) Prediction of the performance of the HSA at steady-state reference 
cycle conditions. Including working fluid, HSHX, and PICS temperatures 

(4) Analysis of +he auxiliary cooling system (ACS) mode of operation, 
including fill gas selection and PICS temperature prediction 

(5) Evaluation of the effect of the BIPS startup transient and prediction 
of low-cycle fatigue (LCF) life 

(6) Analysis of the HSA fueling transient 

(7) Study of the emergency cooling system (ECS) operation to predict the 
PICS temperature excursions 

(8) Resonant frequency and dynamic loading analysis of the HSA and the 
HSHX 

(9) Preparation of a detailed layout drawing of the HSA 

(10) Recommendations for future work for improvement and verification of 
the HSA design. 
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THERMAL ANALYSIS DESCRIPTION ANO DISCUSSION OF RESULTS 


Heat Source Assembly Thermal Model 


The AIRLA general steady-state and transient thermal analyzer computer 
program, H0910, was utilized to perform the required thermal analysis of the 
HSA. The program uses the basic lows of conduction, convection, and radiation 
heat transfer, along with the laws of heat storage and the conservation of 
enerqy to solve probl*ns of heat transport and storage. It employs the elec- 
trlcal analog technique to convert a thermal system of arbitrary complexity in o 
a network model consisting of resistance and capacitance elements or nodes. The 
steedv-state solutions of the thermal model are calculated by the method of suc- 
cessive substitution (Gauss-Sei del method), while the transient solutions are 
calculated by an explicit finite-difference technique. The detailed description 
of the program and the methods of solution are presented in appendix A. The 
following are some of the important features of the program: 


(1) Calculates thermal resistance and capacitance from the inputs of 
system geometry and physical properties. 

(2) Calculates the convection heat transfer coefficient from the input 
fluid properties and the flow characteristics. 

(3) Reevaluates the fluid properties at specified time intervals based 
on updated fi lm temperature. 


(4) Calculates the fluid stream pressure drop. 

(5) Thermal conductivity and specific heat of The material may be input 
as a function of node temperature. 

(6) Heat sources may be Input as a function of time or node temperatures. 

(7) Fluid stream weight flow or fluid velocity may be input as a 
function of time. 


A detailed thermal rodel of the entire HSA was set up tor use with the 
thermal analyzer canputer program. In setting up the model, each component is 
divided into an appropriate number of elements in each dimension, and each e e- 
ment is denoted by a node number. The geometric properties (length, area etc.) 
and physical properties (density, specific heat, and thermal conductivity) of 
each element, as well as the heat transfer mechanisms (conduction, convection 
and radiation) between the adjacent elements, are Input to the program for ca 
cu I at ion of the resistance-capacitance network values. The program calculates 
the temperature distribution and fluid pressure drop of the hSA and HSHX for 
the Imposed boundary conditions and specific configuration. 


The retailed thermal model was used for many of the analysis tasks during 
the study. Modifications were made, where necessary, to reflect the particular 
requirements of each task. The thermal model was essentially two-dimensional 
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with elements in the axial and radial planes of the HSA. Elements in the third 
dimension, the HSA circumferential plane, were added where appropriate. The 
models will be discussed in detail for each particular task. 

Heat Source Insulation System 

The H$IS is designed to limit heat losses from the HSA. The HSIS consists 
of Multi-Foil insulation— numerous parallel layers of a metal I ic foil inter- 
spaced with small particles of a refractory oxide. The foil thickness is 
approx imatel y 0.013 mm (0.0005 in.), and the spacing is around 0.10 to 0.18 ron 
(0.004 to 0.007 in.). The oxide particles prevent adjacent layers of foil from 
contacting each other. The particles offer a low-conductivity thermal path, 
so the Multi-Foil insulation tends to exhibit the thermal conductivity of the 
gas between the layers. The thin metal foils are c' very low resistance. The 
foil layers are wrapped on a mandrel outside of the HSHX. 

When the space between the layers is evacuated, the conductivity path is 
only through the oxide particles that contact successive foil layers. Contact 
resistances between the particles and the foil are quite high. Radiative heat 
transfer between layers becomes significant. At high temperatures, the amount 
of heat transferred by radiation is much greater than that transferred by 
conduction; the Multi-Foil insulation thus resembles a series of radiation 
srx elds. At lower temperatures, both radiation and conduction are important. 

In the thermal analysis, the heat transfer perpendicular to the plane of 
the foils is modeled as an effective radiation using the heat flux correlation 
developed cmpirica! ly by the Thermo Electron Corporation (TECO), the manufacturer 
of the foils (BIPS Preliminary Design Review, AIRPHX, July 1976, p. 4.2-10). 


1.06 x 10 -’ 2 (T 4 - T s|nk <) 
0-778 N + (1 .11 x 10‘ 2 ) N 2 


( 1 ) 


where <p is the heat flux perpendicular to the foils, W/cm 2 
T is the heat source temperature, K 
t S!NK is +he effective sink temperature seen by the last cold foil, K 
N is the number of foil layers 

Due to the low heat fluxes and large number of foil layers, the heat 
source temperature in eq. (I) can be identified with the temperature of the 
hottest foil and the heat sink temperature can be identified with the cold- 
est foil * ith little loss of accuracy. Thus for a given number of foils: 

* * ° (< - T c) 
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( 2 ) 



where o Is the Stefan-Boltzmann constant 
Th Is the hottest toll temperature 
T c Is the coldest foil temperature 

F | 2 Is the effective radiation exchange factor for the entire 
foil packet 

The use of eq. (2) greatly simplifies the thermal modeling because the 
individual foil layers need not be divided Into nodal elements. Heat transfer 
parallel to the plane of the foils was not modeled in the thermal analysis. 

The Multi-Foil insulation for the HSA is in three separate assemblies— 
a cylindrical section and two planar end sections. The interface between the 
sections is a potential heat leak region. In most corner designs, a series of 
gaps exists at the intersection of the foil assemblies. These gaps act as 
cavities for radiation heat leak. TECO edge loss data (BlPS Preliminary Design 
Review, AIRPHX, July, 1976, p. 4.2-11) were utilized in the thermal model to 
account for the heat leakage. The data were treated In a manner similar to the 
perpendicular heat flux. Edge losses also were included for the insulat'on 
discontinuity In the region of the BlPS ducting. The total heat leak from the 
HSA can be roughly broken down as follows: 50 percent from the ends, 30 percent 

from the circumferential surface, and 20 percent from the insulation edges. 

Configuration Study and Candidate Selection 

The configuration study effort analyzed candidate HSHX designs In suffi- 
cient detail to allow comparative evaluation and selection of the configuration 
of greatest merit. Three basic designs were analyzed for each of two reference 
cycles representing a high and a low HSHX working fluid outlet temperature. The 
fluid outlet temperature can be identified with the turbine inlet temperature 
(TIT); actually there is a difference of a few degrees due to duct heat losses. 

The three HSHX configurations chosen for detailed analysis were the 
single-tube helical coil, the helical-flow tube bank, and the axial-flow tube 
bank. Those are shown in fig. 1. Other configurations were investigated, but 
not studied in detail due to what were considered to be inherent disadvantages. 
Some of these alternate concepts are shown In fig. 2. The multitube box mani- 
fold configuration has a large "lost area" region, where the IHS would be 
radiating heat directly to the HSIS. This would tend to increase the HSA heat 
losses and also Increase the HSHX maximum wall temperature. The multitube tan- 
gential start configuration exhibits more effective use of the heat transfer 
surface but would seem to require a quite complex fabrication t«x:hnique with an 
attendant reduction in reliability. The double-row configuration has a quite 
ineffective use of the heat transfer area; the tubes facing the IHS will tend 
to run considerably hotter than the tubes facing the HSIS. 

Mater ial selection . — The proposed material of construction for the high-TIT 
design wa> the refractory columbium alloy C-103. The low-TIT design specified a 
superalloy material. Five superalloys were considered as cand i dates— Hastel loy 
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SINGLE-TUBE HELICAL COIL 



F i ojure 


1. --Candidate HSHX Con f i (jurat ions 
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Figure 2. --Alternate HSHX Concepts. 
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long-term creep data for Hastelloy X. The material has been used successfully 
for more than two decades in a variety of elevated temperature applications 
requiring high strength. Waspaloy represents a potentially attractive alter- 
native exhibiting superior structural properties in the temperature range of 
Interest. However, long-term creep data are sparse and extend for only several 
hundred hours. The lack of a satisfactory data base precludes further consid- 
eration of Waspaloy as an HSHX material at this time. 

The HSHX will be exposed to a space vacuum during operation. Under high 
temperature and vacuum conditions, evaporation of the more volatile components 
of Hastelloy X, such as chromium, may cause subsurface porosity to develop after 
a period of time. Changes in surface composition can cause formation of com- 
pounds and phases that affect surface ductility and consequently fatigue resist- 
ance. Based on current knowledge of the evaporation behavior of Hastelloy X, a 
maximum operating temperature of 788°C (1450°F) has been specified. No maximum 
operating temperature is specified by the C- 1 03 designs. 

Design conditions . — The reference cycle conditions used for the configura- 
tion study are presented in table 3. As previously mentioned, separate TIT 
levels were used for the two HSHX materials. It should be noted that the cycle 
conditions are for a single HSA in a two-HSA BIPS configuration while the pres- 
sure containment requirement Is for a three-HSA BIPS configuration. The working 
fluid is the specified BIPS Xe-He mixture with a molecular weight of 83.8. The 
overall axial length of the HSHX is fixed by the BIPF duct-to-duct centerline 
distance of 40.34 cm (15.88 in.). The assumed surface emissivities represent 
untreated C-103 and roughened and ox?dl 2 ed Hastelloy X. 

TABLE 3 

DESIGN POINT FOR HSHX CONFIGURATION STUDY 



C-103 

Haste 1 loy X 

F low rate, g/s ( 1 b/s) 

44.18 (0.0974) 

53.02 (0.1169) 

Inlet temperature, °C (°F) 

677 (1250) 

— 

Minimum outlet temperature, *C (°F) 

871 (1600) 

704 (1300) 

Maximum metal temperature, *C ( 0 F ) 

— 

^88 (1450) 

Inlet pressure, MPa (psia) 

0.3796 (55.05) 

0.434 (63.0) 

Pressure containment requirement, 
MPa (psia) 

0.793 (115) 

0.793 (115) 

Pressure drop goal (0,3 percent 
AP/P) , kPa (psi) 

1.14 (0.165) 

1.30 (0.189) 

HSHX surface emlssivity 

0.4 

0.8 
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Different approaches were followed In HSHX design for the two materials 
of construction. For the C-103, a constant fluid Inlet temperature was assumed 
with the goal of minimizing the HSHX wall temperature. For the Hastel toy X, a 
maximum wall temperature was assumed with the goal of maximizing the HSHX fluid 
outlet temperature. It should be noted that the fluid temperature rise in the 
HSHX (and hence, heat transfer rate) Is largely dependent on the insulation 
effectiveness and is relatively insensitive to the heat exchanger design. For 
the configuration study, the HSIS was assumed to consist of 110 layers of Multi- 
Foil insulation. 

Axial-flow tube bank .— In the axial-flow HSHX, a series of straight tubes 
is located around the circumference of the IHS. Toroidal type manifolds are 
required, as shown in fig. I. The inner toroidal surface has a few mm clearance 
around the IHS. Sufficient land is required between tubes so that they can be 
attached to the manifolds and maintain structural Integrity. The tube separation 
results in a direct radiation path between the IHS and the HSIS. 

The manifolds are the cause of some of the disadvantages and problems 
associated with this design. It is necessary to carefully size the manifolds 
to minimize flow maldistribution. Sizing was performed using the AIRLA manifold 
design computer program X1000, described in appendix B. Flow distribution con- 
siderations are discussed later in this section. The manifolds tend to use up 
a significant portion of the allowable working fluid pressure drop, while pro- 
viding little additional heat transfer. In addition to wasting pressure drop, 
this situation intensifies flow distribution problems since maldistribution 
is closely related to the ratio of manifold pressure losses to heat exchanqer 
core pressure drop. 

The poor heat transfer conductance in the manifolds results in high 
manifold metal temperatures, with the outlet manifold becoming the HSHX hot 
spot. For the Hastel loy X design, the outlet manifold was assumed to be insu- 
lated, bringing the manifold temperature below the maximum tube wall tempera- 
ture. The insulation was assumed to be 10 layers of Multi-Foil shielding. 

The thermal model constructed for the axial-flow tube bank is shown in 
fig. 3. Only one heat exchanger tube Is shown; the wall temperatures and work- 
ing fluid t emperatures and flows are assumed identical for all tubes. The tube 
shown is n radians from the duct location, at the region of lowest velocity in 
the manifold. Each tube is divided Into two radial sections~an Inner section 
that faces the IHS, and an outer section that faces the insulation. Fluid nodes 
are shown In the center of the tube and In the manifolds. The outlet manifold 
Insulation is not shown, but is Incorporated In the model for the Hastel loy X 

fjrj r I «« ' 

U«a)l i • 


Since the PICS temperatures were not investigated In the configuration 
study task, a simplified model of the IHS was used. The model consisted only 
of Internal and surface nodes. The two IHS end nodes represent the laminated 
end crush-ups and the graphite support plates. The support plates are not 
physical I \ part of the IHS. 
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HOPES COMPONENT 

iT" HEAT SOURCE INTERNAL 

9-16 GRAPHITE AERO-SHELL 

17-28 "ASTRO -QUARTZ" INSULATION 

29-36 ZIRCON IA ENO SPACERS 

J7-AO INCONEL X-7S0 PRCLQAO NUT 

41-62 TITANIUM SUPPORT SPIDERS 

43-66 TITANIUM SUPPORT FLANGE 

65-56 LOCKALLOY HOUSING 

57-80 MULTI-FOIL INSULATION <110 LAYERS) 

81-96 INNER HSHX SURFACE 

95-108 OUTER HSHX SURFACE 

I09-H7 INLET MANIFOLD 

1IJ-II6 OUTLET MANIFOLO 

200-213 Xr-Hc WORKING FLUID 


« I 4J 



"* 85 . 


12 - 87 . 


'> * 89. 


'S • *K 


NK. 






Figure 3.— Multitube HSA Thermal Model. 
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Heat Is transferred by radiation from the IHS to the HSHX, from the HSHX 
to the HSIS, and from the IHS to the HSIS through the gaps between the tubes. 
Radiation exchange factors were calculated for the different connections. 

Heat transfer In the HSHX Is characterized by conduction In the tube and 
manifold walls, Internal radiation from the Inner tube node to the outer tube 
node, and forced convection between the working fluid and the walls. In the 
manifolds, the fluid velocity decreases with circumferential distance from the 
entrance (inlet manifold) or exit (outlet manifold) port. Thus, the heat trans- 
fer coefficient will also decrease. The thermal model considers the location 
tt radians from the ports, where the velocity Is lowest and the heat transfer 
poorest. Because the fluid is almost stagnant, conduction should be the pri- 
mary mode of heat transfer. This conservative approach yields the minimum 
fluid temperature Increase and the maximum manifold wall temperature. 

Pressure drop in the manifold region inlcudes losses for transitions 
between the connecting ducts and the manifolds, flow frictional losses in the 
manifolds, and turning losses between the manifolds and the heat exchanger 
tubes. These losses are calculated with the aid of the AIRLA manifold design 
computer program, XI000 (appendix B). 

The heat transfer and pressure drop for the tubes are calculated by the 
thermal analyzer program from Input values of the tube Colburn modulus and 
Fanning friction factor, respectively. In order to enhance the heat transfer 
characteristics of the axial tubes, ring dimples are assumed in the model. 

AIRLA has fabricated ring-dimpled tubes of various diameters and materials. 

The dimples are formed in the tube wall by a rolling process that changes the 
tube wall structure only slightly. Dimple geometry can be varied in both spac- 
ing and depth over wide ranges. A sketch of a ring-dimpled tube is shown In 
fig. 4. As a point of comparison, the axial tube bank ring dimple design has 
a Reynolds number of around 2000. At this Reynolds number, the dimpled tube 
Colburn nod u I us is about twice that of a smooth tube. Ring dimpling increases 
the pressure drop by a factor of about 2.7 over that for a smooth tube. 



Figure 4.— Ring-Dimp led Tube Geometry. 
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Tdblos 4 and 5 summarize the seml-optlmlzed axial-flow tuba bank designs, 
tba C-103 and Hastalloy X conf I jurat l on* a ra similar, axe apt for the out I at 
f manifold Insulation for tha Hastatloy X case. Tha staady-stata MSA temperature 
prof lias art shown In figs. 5 and 6 for tha C-103 and Hastalloy X daslgns, 
respectively. Tha Innar and outar surfaca tamparaturas along with tha working 
fluid tamparaturas a ra shown for tha hast exchanger. Tha manifold matal tam- 
paraturas a ra shown; however, tha manifold Insulation tamparaturas for tha 
• Msstat loy X configuration a ra not shown. Tha hast axchangars consist of 83 
rlng-dlmpled tubas, 6.33 mm (0.23 In.) In dlamatar. Tha tuba and manifold wall 
thickness raqulramants wars datarmlnad by a prassura contalnmant analysis, which 
Is discussed in tha structural analysis section. 

Both daslgns meat tha fluid outlet tamparatura raqulramants of 8?t*C 
(1600*F) for tha C-103 HSHX and 704*C (1300*F> for tha Hastalloy X HSHX. The 
prassura drop goals are achieved with tha manifolds using about 30 percent of 
tha al lowabta. 

Hellcal-flos tuba bank .— In tha hal leal-flow tuba bank, a series of tubas 
bant Into a helical coll so that each tuba makes exactly one complete turn 
around the IMS. Tha toroidal manifolds require spe^e b e twee n tha tubas In tha 
mm 1 fold region for tuba attachment. 

If tha tubas are colled Into a helix of a particular angle, however, the 
. tubas will be In contact over most of tha heat exchanger length as can be seen 
In fig. 1. Tha required angle Is a function of tha number of tubas usad and 
tha tuba dlamatar-to-col I dlamatar ratio. With tha tubas In contact, essen- 
tially all of tha radiation emitted from tha IMS In tha radial direction Is 
absorbed by tha HSHX. 

Tha tuba curvature enhances tha heat transfer and Increases tha prassura 
drop. Because tha tubas are colled, tha actual flow length is considerably 
greater than for tha straight-tuba configuration. This also tends to Increase 
both tha heat transfer conductance and tha pressure drop. Thus, the tube diam- 
eters are somewhat larger than for tha straight tuba In order to meet the 
prassura drop requirement. 

Tha manifold considerations are similar to those for the straight tube 
bank. To attain flow uniformity by balancing tha manifold areas, it is neces- 
sary for each tuba to have tha same relative position In both the Inlet and 
outlet manifolds. This requires on integral number of turns. A single turn 
was chosen, due to the limited pressure drop allowance. 

The thermal model for the helical-flow tube bank Is similar to that for 
the axial-flow tube bank (see fig., 3). Although HSHX symmetry is not precise 
In the HSA circumferential direction, the temperature variation In this direc- 
tion dhould be quite smell compared to the axial temperature variation. Thus, 
the two-dimensional HSHX model was deemed suitable for the helical tube bank. 

The actual curved tube Is represented by a straight tube In the thermal model; 
the actual tube flow length is used. 
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TABLE 4 

axial-flow tube BANK DESIGN 

Number of tubes 
Tube OD, mm (in.) 

Oimple spacing, mm (In.) 

Dimple depth, mm (in.) 

,ln ® Tube diameter, cm (In.) 

Tube length, cm (In.) 

Inlet manifold ID, nn (in.) 

n n i?L manffold " aM Thickness, tm (In.) 

Outlet manifold ID, inn (In.) 

Outlet manifold wall thickness, mm (In.) 


83 
6.35 
9.53 
0.38 
8.89 

23.5 

36.6 
27.8 
0.33 
41.1 
1.40 


(0.25) 

(0.375) 

(0.015) 

(0.35) 

(9.26) 

(14.40) 

(1.095) 

(0.013) 

(1.62) 

(0.055) 


■ ndLfc 5 


AXIAL-FLOW TUBE BANK PERFORMANCE 


Tube wall thickness, mm (In.) 

Weight* (Including manifolds), kg (lb) 
Fluid Inlet temperature, °C (*F) 

Fluid outlet temperature, *C (*F) 

Maximum tube wall temperature, *C CF) 
Maximum manifold temperature, *C (*F) 

Tube pressure drop, kPa (psl) 

Total pressure drop, kP a (psl) 

circumferential tub. gradient, 

HSA heat loss, percent 
Maximum INS surface temperature, *C CF) 

•Includes a wrap-up weight of 0.45 kg (1 lb) 


C-103 

Haste 1 loy X 

0.254 (0.010) 

0.254 (0.010) 

3.2 (7.1) 

3.0 (6.6) 

677 (1250) 

577 (1070) 

884 (1623) 

752 (1386) 

931 (1708) 

787 (J449) 

941 (1725) 

778 (1443) 

0.538 (0.078) 

0.565 (0.082) 

1.14 (0.165) 

1.23 (0.179) 

10 (18) 

20 (36) 

5.2 

3.7 

959 (1759) 

847 (1557) 
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TEMPERATURES IN °C 


TEMPERATURES IN °R 



Mi- 



*> 



Mi 



•-2M27 


NOTES s V « 44.18 g/SEC (0.0974 LB/SEC) 
P - 0.3796 MPt ( 55.05 PSIA) 

SINK TEMPERATURE - 49°C (580°R) 





Figure 5.-- 0103 Axial-Flow Tube Bank HSA Temperature Profile. 
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pr««r. dl^p ^ * Thl*he r I « ^* C h r V Sr U r" °" hMt ,r< "’ s, «‘ and 
°* f 0 *' 1 3500. This Is well .itMn "1 !s»?„^"I, d * Sl9n has a R «y"°'«s 
tub« curvature Increases the transition R«#n«?rtl ! f ow ro 9 |me since the 
-ra«. !> 9 lv. th. transition R^l'dT^^ ~ 


R *t ■ 2100 Cl + 12 (d/D)0*5j 

curvature* T, ^sl^on^R^'ds^nu^ “‘around ^ d “ , S" 

Varrler <?£?!“)* ^used? ^ ,0minar reg,me * the correlation of Kubair and 


(3) 


(4) 


f = 0.7716 Cexp (3.553 d/D)3R e -0«5 
where f Is the Fanning friction factor. 

a Irict Iw^i^oMuIul?" V*’ str^t^ ■T"?*? °' 350 °- «• < 4) Predicts 
number Is around 0. 007.* * straight- tube friction factor at this Reynolds 

tor fuMy^velop^^i^ar nJTn'uIid^Thl'f l° # Naka >' a ™ (ref. ?) 

eloped close to the inlet becni« n +h * The f ow ,s * ikel V to be fully dev- 
and hydrodynamic entry length (ret. J): Ur ”l U cor?^t,o„ r ^ UC ' n9 the ,h! - r ™ 1 


Nu/Nu 


0.1978 0.5 

7 On 



•mi 


( 5 ) 


( 6 ) 


where Nu is the Nusselt number 


h.a% fiuVtt ,r ?u?rr^iop« e !L?“^; M wu ‘' ,o 4 - j64 ,or uni,or '" 

Pr is the Prandtl number 

— layer thinnesses. 

On Is the Dean number given by 

On = Re < d/D ) 0 . 5 

( 7 ) 
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The Colburn modulus, j, Is related to the Nusselt number by: 

j ® NuRe“^ Pr“^ 


As a point of comparison, at a Reynolds number of 3500, eqs. (5) through 
(8) predict a Colburn modulus of 0.00347. The straight-tube Colburn modulus 
at this Reynolds number is around 0.0023. 

Tables 6 and 7 summarize the semi-opt imized helical-flow tube bank designs. 
The C- 103 and Hantelloy X configurations are similar, except for the outlet 
manifold insulation for the Hastel Joy X case. The steady-state HSA temperature 
profiles are shown in figs. 7 and e for the C- 103 and Hastel I oy X designs, 
respectively. A single tube of the heat exchanger is shown. For convenience, 
the tube is viewed straight, although the actual tube is curved and of much 
greater length than shown. The inner and outer surface temperatures along with 
the working fluid temperatures are shown for the heat exchanger. The manifold 
metal temperatures are shown; however, the manifold insulation temperatures for 
the Hastelloy X configuration are not shown. The heat exchanger designs consist 
of 39 tubes, 7.94 mm (0.3125 In.) in diameter, each making one turn around the 
IHS* A conceptua I design drawing is shown in fig. 9. The tube and manifold 
wall thickness requirements were determined by a pressure containment analysis, 
which is discussed in the structural analysis section. 

Both designs meet the fluid outlet temperature requirements of 871°C 
( 1600°F) for the C-103 HSHX and 704 8 C (1300°F) for the Hastelloy X HSHX. The 
pressure drop is somewhat greater than the design goal. The manifolds use 
close to 50 percent of the pressure drop. 

SJ.hgl e-tube h elical coil .— In this design, a single tube is wrapped around 
the IHS in the shape of a tightly wend helix,. The separation between successive 
turns of the coil is very small so 1r,*t very little radiation is transferred 
directly from the IHS to the HSIS. The single tube has essentially no flow 
distribution problem and negl igible manifold pressure drop. Thus, all of the 
pressure drop can be expended in the heat exchanger proper to yield maximum 
feat transfer conductance. The relatively severe tube curvature enhances the 
heat transfer coefficient and increases the pressure drop. See fig. 1 for a 
sketch of this configuration. 

The HSHX manifolds are essentially extensions of the heat exchanger itself 
with high heat transfer conductances. Thus, the manifolds need not be insulated 
to meet the maximum wall temperature requirement for the Hastelloy X design. 


For the single-helical-tube heat exchanger there is a significant variation 
in +emperature profile in the circumferentia I direction. A modified three- 
dimensional model is used for thermal analysis. The IHS surface and the heat 
exchanger are divided into four circumferential nodes, in addition to the radial 
and axial nodes similar to the multitube model. The remainder of the HSA model 
is essentially two-dimensional and Identical to the multitube configuration. 
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TABLE 6 

HELICAL-FLOW TUBE BANK DESIGN 


Number of tubes 
Number of turns 
Tube 00, mm ( in. ) 

Center 1 1 ne col I diameter, cm (In.) 

Total tube length, cm (In.) 

Helix angle, radians 
Inlet manifold ID, m (In.) 

Inlet manifold wall thickness, mm (in.) 
Outlet manifold ID, nvn (In.) 

Outlet manifold wall thickness, mm (In.) 


39 

1 

7.94 (0.3125) 
23.5 (9.26) 

83.8 (33.0) 
0.433 

27.8 (1.095) 
0.36 (0.014) 
41.1 (1.62) 
1.40 (0.055) 


TABLE 7 


HELICAL-FLOW TUBE BANK PERFORMANCE 


Tube wall thickness, mm (in.) 

Weight* (including manifolds), kg (lb) 
Fluid Inlet temperature, *C (*F) 

Fluid outlet temperature, *C (*F) 
Maximum tube wall temperature, *C (*F) 

| Max I mum manifold temperature, *C (*F) 
Tube pressure drop, kPa (psi) 

Total pressure drop, kPa (ps!) 

Maximum circumferential tube gradient, 

HSA heat loss, percent 

Maximum I HS surface temperature, *C (*F) 


C-103 

Hastel loy X 

0.279 (0.011) 

0.254 (0.010) 

3.9 (8.5) 

3.4 (7.5) 

677 (1250) 

542 ( 1 0C 7 ) 

884 (1623) 

718 (1324) 

951 (1744) 

788 (1450) 

958 (1757) 

757 (1394) 

0.772 (0.112) 

0.745 (0.108) 

1.37 (0.199) 

1.39 (0.201) 

25 (45) 

43 (78) 

5.3 

975 (1787) 

3.5 

837 (1538) 


•Includes a wrap-up weight of 0.45 kg (1 lb) 
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TEMPERATURES IN R 


W «= 44.18 g/SEC (0.097 1 * LB/SEC) 
P * 0.3796 MPa (55.05 PS IA) 

SINK TEMPERATURE - 4g°C (580°R) 


S 70049 


Figure 7 . — —C— 103 Helical-Flow Tube Bank HSA Temperature Profile. 
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TABLE 8 


SINGLE-TUBE HELICAL COIL DESIGN 


Number of turns 

10 

Tube 00, mm (In.) 

43.18 (1.70) 

Centerline coll diameter, cm (In.) 

23.9 (9.41) 

Total tube length, cm (In.) 

762 (300) 


TABLE 9 

SINGLE-TUBE HELICAL COIL PERFORMANCE 



C-103 

Haste 1 loy X 

Tube wall thickness, wn (in.) 

1.09 (0.043) 

1.35 (0.053) 

Tube weight, kg ( lb) 

9.75 (21.5) 

11.1 (24.4) 

Fluid inlet temperature, *C (*F) 

677 (1250) 

578 (1073) 

Fluid outlet temperature, # C ( F) 

886 (1626) 

754 (1390) 

Pressure drop, kPa (psi) 

1.13 (0.164) 

1.21 (0.175) 

Maximum metal temperature, °C (*F) 

916 (1680) 

786 (1447) 

Maximum circumferential tube gradient, *C (*F) 

52 (93) 

56 (101) 

HSA heat loss, percent 

4.6 

3.4 

Maximum IHS surface temperature, *C (*F) 

960 (1760) 

824 (1515) 







TEMPERATURE IN °R 


TEMPERATURES IN U C 


W - 44.18 g/SEC (0.097** LB/SEC) 
P - 0.3796 MPa (55.05 PS IA) 

Tin - 677°C (1710°R) 

SINK TEMPERATURE - 49°C (580°R) 


Figure 11.--C-103 Single-Tube Helical Coll MSA Temperature Profile 
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w - 53.02 g/SEC (0.1169 LB/SEC) 
P = 0.434 MPa (63-0 PS IA) 

Tin ■ 579°C (1533°R) 

SINK TEMPERATURE - 49°C (580°R) 


Figure 12.— Haste) toy X Single-Tube Helical Coil HSA Temperature Profile 


















Flow distribution .— The attainment of uniform flow distribution is an 
important consideration for the multitube configurations. Flow maldistribution 
will result in higher surface temperatures In the tubes with low mass flow. In 
addition to the direct effect of the higher metal temperature, the the-mal gradi- 
ents will be increased, thereby influencing the structural integrity of the unit. 
Also, flow maldistribution may cause an increase in the overall pressure drop. 

The flow distribution is controlled by the flow areas of the manifolds. 

This is illustrated in the following discussion. Consider fig. 13, which 
shows typical heat exchanger manifold configurations. Fig. 13a shows that 
all elements of the flow stream experience a common pressure at the heat 
exchanger inlet, Pj, and a common pressure at the heat exchanger outlet, Po. 
Consequently, the pressure loss through all flow streams must be identical, 
i.e., P] - If two flow paths have different resistance characteristics 

the flow Is apportioned as necessary to balance the overall pressure drop. 

Thus, in fig. 13a, the remote path A will receive a smaller portion of the 
total flow than the more direct path B. 

Pressure losses fall into two categories: (1) frictional (or viscous 

shear) and (2) kinetic energy. In practical heat exchanger design, frictional 
losses are usually predominant In the core (or heat transfer matrix) and 
kinetic energy losses predominate in the manifolds; however, all losses must 
be considered for optimum design. As a general design guideline, reasonably 
uniform flow distribution can be achieved by designing manifolds such that 
their loss character i st ics are low, relative to the frictional loss character- 
istics of the core. The disadvantage of this is that for some configurations, 
the physical size of the manifolds must be larger than would be the case if 
detai led analysis were applied. Consider, for example, cases in which the 
oriertation of supply and discharge lines must be perpendicular, rather than 
parallel, to the direction of the flow through the heat exchanger core. Typical 
■samples are shown in fig. 13b. In both of these cases, the cross-sectional 
manifold area available for distribution of the flow is proportional to the 
small dimension Y rather than, as in the case of parallel supply and discharge, 
the large dimension X. Since both frictional and kinetic energy losses are 
proportional to the square of the flow stream velocity, this has a profound 
effect on the manifold losses. Unless core frictional losses are very high, 
the configuration shown In fig. 13b is inherently inferior to that shown in* 
fig. 13a. 

In comparing the configurations of fig. 13b, the initial impression is 
that di st»- ibution of the Z-f low configuration is superior to that of the U-flow 
because of the latter's disparity with respect to the ratio of extreme flow 
path lengths. In Z-f low, all flow streams have the same path length. In fact 
however, U-flow is superior to Z-f low. This Is due to a factor heretofore 
neglected, i.e., the conversion of velocity pressure to static pressure. It 
is observed in fig. 13 that the inlet velocity pressure acts in the direction 
of core flow; it Is therefore available to assist the circulation directly. 
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For constant-area, un i form— f low manifolds, the inlet manifold velocity 
pressure is much larger than the cutlet manifold velocity pressure. Thus, if 
the density change is not too great, the inlet manifold flow area will be’ 
smaller than the outlet manifold area. With this in mind, the design tech- 
nique involves assigning an area for the outlet manifold (such that the 
envelope requirements are met) and finding the best inlet manifold area. 

The AIWLA manifold design computer program, X1000, was used for manifold 
sizing. The calculations include the effects caused by frictional losses in 
the manifolds. Because the toroidal manifolds are curved, friction factors 
for flow through curved tubes are used. Due to the decreasing flow velocity, 
the flow is turbulent near the mainifold open end and laminar near the closed 
end. In the turbulent region, Ito's correlation, eq. (92, is used. In the 
laminar region. White's empirical formula as given by I to (ref. 5) is used. 

The relatively severe manifold curvature is beyond the range of the correlation 
of Kubair and Varrier (ref. 2). 

The overall pressure drop Is determined by summing the heat exchanger 
core losses, the manifold turning and frictional losses predicted usingcomp- 
puter program X100C, and the transition losses between the manifolds and the 
connecting ducts. The latter are quite significant due to the high losses 
associated with diverging and converging branched flow (ref. 6). Computer 
program X1000 is described in appendix B. 

Candidate comparison and selection .— A comparision of the candidate 
designs is presented in table 10. It can be seen that the axial-flow tube bank 
is clearly the better of the two multitube designs. The weight and thermal per- 
formance are superior to the helical-flow tube bank. The HSIS is so effective 
that the direct impingement of radiation from th R IHS on the HSIS in the axial- 
flow design does not substantially increase the thermal losses. The heat is 
essentially re-radiated to the back side of the HSHX by the insulation. For the 
C-103 design, the axial tube bank heat leak is actually lower than the helical 
tube bank heat leak. This is a result of the lower HSHX wall temperature due 
to superior heat transfer for the ring-dimpled axial tube bank. 

A comparision of the axial-flow tube bank with the single-tube helical coil 
indicates substantially similar thermal performance. The axial-flow maximum 
metal temperaf e for the C-103 design exists on the outlet manifold, which is 
not insulated. Use of insulation as in the Hastel loy X design would lower the 
maximum metal temperature. 

It can be seen that the axial-flow tube bank offers a substantial weight 
savings ever the single-tube helical coil. Nevertheless, the selected design 
is the single-tube helical coil. The inherent simplicity and potential struc- 
tural integrity and reliability of the single-tube design are considered to be 
of greater importance than the minization of weight. The multitube design 
requires over one hundred tube-to-man i fold welds representing a large number of 
potential problem areas. The single-tube design has other advantages including 
no flow distribution problems, no requirement for manifold insulation, and neq- 
I i gi bl e man i fol d pressure drop. s 
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table 10 


hshx weights Af.o perforwamce 


height, kg ( , b) 

T'uiO outlet temperature, X CF> 

HSA heat loss, percwf 

t'aximue metal temperature. *c f F , 
Hastef Iny y 
"eight, kg (i b ) 

Fluid outlet tenperature, °C (°F> 

** Teat loss, percert 

v a*fmum metal temperature, ’C CF) 


S i nc I e—tu be 
h el ical coil 


9.75 ( 21 . 5 ) 
686 (1626) 


916 (1680) 


H.l (24.4) 
754 (1390; 


786 (1447) 


Axial -f low 
tube bank 


3.2 (7.1) 
884 ( 1623 ) 


941 (1725) 


3.0 (6.6) 
752 (1386) 


787 ( 1449 ) 


Hel ical-ficw 
tube bank 


3.9 (8.5) 
884 (1623) 


958 (1757) 


3-4 (7.5) 

7 f 8 (1324) 


788 (1450) 


, . _ ax la I -t low tub® bonk to bo a viable alternative. 

„ ,.*£? r TWit 

...Hi tub. bank shou'ld considered. 

Hastel loy X I. the 

emphasis on MSA rel 1 1 ty , Ai eor . ^ ^ difficult to limit the P'^$ 

Prol Imlnary analysis Indicate* ]\ a , tor tht , h | q h fluid temperatures In 

TOW “cr.-W irmirr^ul-ents can bo' met at the lower operat- 

tenure o. .be Hastel loy X design. 


Meat Source Heat Exchanger Optimization 


... iichy coi.xrt«fl as the slnql e-tube helical 
With the configuration ot 'e - - th€ hM t exchanqer »os under- 

coll Hastel loy X design, the task ot of m f T ^ n.sed on Into, motion 

taken. The design reference cyel. wndltlw. **. Melons, which replace 

from the RIPS program of o table II This reference cycle was used 

those given In table 5. are given In table K the tw-HSA 

tor the remainder ot the study. The * ™en, I s tor the three-HS» 

configuration while the pressure containment regulreme 

design . 

TABLE II 

DFSIGN REFERENCE CYCLE conditions 


F low rate, g/s ( Ib/s) 

Minimum outlet temperature, “C ( # l) 

Maximum metal temperature, *C Cn 
• n 1 1 1 pressure, MPa (psia) 

Pressure containment requirement, MPa (psia' 

Pressure drop goal , percent At /l 

HSHX surface emlsslvlty 

Working fluid temperature rise, *L ( II 


si .778 (0.H415) 
704 (1300) 

788 ( 1 4 so ) 
0.4?67 (M.8«) 
0.703 ( 10 ?) 

< 1 
0.8 

17 ? ( 310 ) 


The limitation considered HSHX ^^"w^X Ve^ge . 

tormenre. Since the HSA hen losses nr ... ntl vely ' Nd , ture. The 

rad^.'Jrtlrh^r^Telrht/^ur/drop. and HUN out, .t temperature 


3b 


1 ! I 


i -l. 


i > —ii 


were acquired from the BIPS program office at AIRPHX (J. Lloyd, BIPS-GDS 
Coordination Memo A0563, November 23, 1977) and are given below: 

AW * C x A(AP/P) (11) 

AW ■ 0 x AT (12) 

where AW Is the change In equivalent weight 

A(AP/P) Is the change In percent pressure drop 
AT Is the change In fluid outlet temperature 
C * -9.07 for AW In kg and -20 for AW In ID 

D * 0.408 for AW in kg and AT in *C and 0.5 for AW In lb and AT In *F 


Eq. (11) Indicates that AP/P could be Increased by 0.2 percent if 1.8 kg (4 lb) 
could be saved In weight. Eq. (1?) Indicates that the fluid outlet temperature 

could be lowered by 0.55°C (1*F) if 0.22 kq (0.5 lb) could be saved In weight. 

The value for D is at a fluid outlet temperature of 760*C ( 1 4 00 " F ) , which will 
be seen to be the approximate value attainable for the optimised design. 

The ground rule used for the HSHX optimization was an integral number of 

turns in the helical coil with on overall linear dimension of 40.34 cm (15.88 

in.), the distance between the centerline of the connecting BIPS ducts. The 
tube outside diameter is equal to the duct centerline separation divided by 
the number of turns— with a 0.51- to 0.7t>-mm (0.02- to 0.03-in.) allowance for 
clearance between successive turns. Tho tube wall thickness requirements were 
determined by a pressure containment analysis, which is discussed in the 
structural analysis section. Thus, the only Independent variable in the HSHX 
optimization is the number of turns. 

HSHX designs are presented In table 1? tor 9, 10, and 11 turns. Colls with 
more than 11 turns had pressure drops exceeding the I percent AP/P restriction, 
while coils with fewer than 9 turns required excessive wall thickness. The 
nominal wall thickness was taken as 5 percent qreater than the minimum allowable 
to account for manufacturing tolerances. A constant clearance was maintained 
between the IHS surface and the inside of the coll; thus the smaller tubes have 
a smaller coil diameter. The HSIS and HSA outer case locations were not changed. 

Heat transfer and pressure drop calculations were performed for the designs 
given In table 12 using the HSA thermal model modified for the specific number 
of turns. HSHX weights were calculated taking Into account the required thick- 
walled ellows (see the structural analysis section). The rosults of the analyses 
art presented In table 13. 
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table j? 

hshx designs for optimization STUDIES 


Number 

of 

turns 


Tube OD 


cm 

in . 

4.42 

1.74 

3.96 

1.56 

3.61 

1.42 


Nominal wai I 
th ickness 


1-42 0.056 
^ - 30 0.051 
1-19 0.047 


Col I 

tJ inmeter 
cm in . 
24.8 9.76 

24.3 9.58 

-4.0 9.44 


Tube 
I enqth 


cm 

in. 

701 

276 

764 

301 

829 

326 


TABLE 13 

HSHX OPTIMIZATION RESULTS 


Numb er 
of 

turns 


AP/P, 

percent 

0.275 

0.504 

0.855 


Fluid outlet 
temperature 


750.1 1382.1 
756.0 1392.8 
760.5 1400.9 


Actual weight 

fcp lb 

11.3 24.9 

10.0 22.2 

9 -1 20.1 


f C|U i valent weight 

kg 

lb 

11.6 

25.7 

10.1 

22.2 

10.5 

23.1 


It c 
and fluid 
?>hou Id be 
turo of 7. 1 
an equiva 
Table 13 
greater, , 
than the 1 
a 1 . 50 -mm 
*ork in ti 


d ou '"'’ C t^P^ature t ^ c "^ r ,h<1 f eSSUre drop 

5 ent weight was determined usino thp I ^ ' PS in «**• HI) and (12) 
■ndlcates that the ll-turn HSHxVts an 0 " 1 "™ f' >5 'S n 05 ’bo baseline. 

• nd the 9-turn HSHX has an equivalent .T-* 8 2\" e ,'f h ! °- 4 k 9 10.9 lb) 
turn unit. Thus the 10-turn HSH* • +? * Greater 

' <0-0 '-.-".I wall Is Ident , M ed ^he Inti” '’TV «> and 

is study concentrates on this conflglr^il A " su ».equeot 


Steady-State Results 


With the selection of the optimized HSHX design, the steady-state HSA 
thermal analysis was performed. The HSA thermal model was updated to Include 
the latest information on the HSA components. The complete thermal model Is 
shown In fig. 15. Table 14 lists the component materials used h the model. 

The number of foil layers in the HSIS was reduced to 60, reflecting the 
lower temperature superalloy HSHX cycle conditions. A Haste I I oy X inner support 
mandrel for the HSIS was included in the model for both the cylindrical and 
)lanar insulation. 

The number of HSHX radial nodes was increased from 2 to 4 per turn. This 
improvement resulted In 160 HSHX nodes, along with 41 fluid nodes. 

A detailed model of the IHS was prepared to enable accurate prediction 
of the PICS temperatures. The IHS model Is strictly two-dimensional, as the 
circumferential nodes on the surface were found to be unnecessary (note the 
small circumferential temperature variation in fig. 12). Within the IHS, heat 
is transferred primarily by radiation, with solid conduction heat transfer where 
appropriate. The partial pressure of helium released by the radioactive decay 
of the plutonium dioxide fuel is assumed to be too low to offer a significant 
conductive path within the IHS. 

A key assumption that has a significant effect on the predicted PICS tem- 
perature is the emissivity of the PICS surface. According to Mound Laboratory 
(personal communication, E.W. Johnson, December 16, 1977), a reasonable value 
for the suface emissivity of the grit-blasted iridium is 0.28. This value was 
assumed for the study. The thermophysical properties for the various forms 
of graphite used in the IHS were obtained from the General Electric Company 
(R.J. Hemler, BIPS-GDS Coordination Memo GO 161, December 12, 1977). 

Using the complete thermal model, a steady-state HSA thermal map was pro- 
duced based on the reference design conditions given in table 11. Pressure 
drop and heat transfer in the HSHX were predicted using eqs. (9) and (10), 
respectively. The average heat transfer coefficient in the heat exchanger tube 
was 108 W/(m^- K) [J9 Btu/(hr-sq ft-*F)]. A heat sink temperature for the HSA 
of 49°C (120°F) was assumed. 

The thermal map Is presen + ed in fig. 16 for all HSA components except 
the HSHX. The HSHX temperatures are listed in tables 1* and 16 for SI and 
customary engineering units, respectively. The HSHX nodal locations are 
Indicated in fig. 17. Each of the 10 heat exchanger turns is divided into 
four quadrants in the coil circumferential (fluid flow) direction and four 
sections in the tube circumferential direction. The fluid temperatures are 
considered uniform within each quadrant. 
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Figure I5.-HSA Thermal Kodel. 



TABLE 14 


HSA COMPONENTS 


Component 

Post-impact containment shell (PICS) 
Graphite impact shell (GIS) 

Retaining ring 
Aeroshel I 
Ablation sleeve 
Comp I i ance pad 
End crush-up 
Aeroshel I end cap 
Support plate 
Fiber insulation 
End spacer 
Preload nut 

End enclosure (spider) 

Support flange 
End dome 

Insulation mandrel 
HSIS 

Outer shel I 
HSHX 


Mater ia I 

Iridium 

Thorne I -50 

Pyrocarb-508 

Poco AXF-5Q 

Pyrocarb-406 

Thornel-50 

Pyrocarb-508 

Poco AXF-5Q 

Poco AXF-5Q 

Astro-quartz 

Z i rcon i a 

Inconel X-750 

T i tan ium 

Titanium 

Locka I loy 

Haste I loy X 

Multi-Foil insulation 

Locka I loy 

Hastel loy X 



TEMPERATURES IN U C 


W = 51-78 g/SEC (0.1142 LB/ 
P = 0.4267 MPa (61 .89 PSIA) 
SINK TEMPERATURE = 49°C (58 


Figure 16. — S+eady-S+ate 
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4 
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1 

2 
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4 
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2 

3 

4 

1 

2 


Fluid 

Out of sect 


1537 

1550 

1562 

1573 

1584 

1593 

1602 

1610 

1618 

1626 

1635 

1ft4? 

1650 

1656 

1666 

1674 

1681 

1689 

1697 

1704 

1711 

1720 

1728 

1735 

1742 

1750 

1757 

1764 

1771 

1778 

1785 

1792 

1798 

1805 

1812 

1816 

1824 

1832 
1840 
184 7 
1853 


1665 

1670 

1676 

1682 

1687 

169? 

1696 

1701 

1718 

1723 

1727 

1731 

1748 

1752 

1756 

1760 

1775 

1780 

1783 

1787 

1807 

1811 

1815 

1818 

1831 

1835 

1839 

184? 

1855 

1858 

1862 

1865 

1876 

1879 

188? 

1886 

1901 

1905 

1908 

1910 


1605 

1613 

1622 

1630 

163? 

1639 

1645 

1651 

1664 

1670 

1676 

1681 

1695 

1701 

1707 

171? 

1725 

1731 

1736 

1741 

1756 

1762 

1767 

1772 

1784 

1789 

1794 

1799 

1811 

1816 

1821 

1826 

1837 

1841 

1846 

1850 

1868 

1872 

1875 

1878 


1596 

1603 

1611 

1618 

1614 

1620 

1626 

1632 

1644 

1650 

1655 

1661 

1675 

1681 

1686 

1691 

1706 

1711 

1716 

1721 

1737 

174? 

1747 

175? 

1766 

1771 

1776 

1780 

1794 

1799 

1803 

1807 

1823 

1827 

1831 

1835 

1860 

1865 

1869 

1873 


Temperature at circumferent I 


B 


a I position 


D 


1609 

1615 

1622 

1629 

1632 

1639 

1645 

1651 

1664 

1670 

1676 

1681 

1695 

1701 

1707 

171? 

1725 

1731 

1736 

1741 

1756 

1762 

1767 

1772 

1784 

1789 

1794 

1799 

1811 
1816 
18 21 
1825 

1837 

1841 

1846 

1850 

1867 

187? 

1877 

1881 


1853 

















The key results of the steady-state analysis are shown In the first column 
of table 17. The PICS temperatures are somewhat below the specified limits of 
1150* to 1200*C for steady-state operation (see table 1). The variation in PICS 
tenperature Is In the HSA axial direction, as can be seen from fig. 17, and is 
a result of the varying fluid temperature In the HSHX and the heat loss from 
the HSA ends. The maximum ablation sleeve temperature, listed for reference, 

Is well within the capability of the ablation sleeve material, Pyrocarb-406. 

The working fluid outlet temperature of 756 *C (1393*F) easMy meets the minimum 
requirement of 704*C (1300*F). The working fluid temperature rise, HSA heat 
loss, and pressure drop meet the reference cycle requirements (see table 11). 

The maximum HSA outer shell temperature Is well within the capability of the 
shell material, Lockalloy. 


TABLE 17 


HSA STEADY-STATE THERMAL ANALYSIS RESULTS 


HSHX emissivity 

0.8 


0.3 

0.25 

P ICS temperature 
range, *C 

1111 to 
1165 

1117 to 
1170 

1138 to 
1188 

1145 to 
1195 

Maximum ablation 

sleeve temperature, 
*C (°F) 

821 (1510) 

829 (1525) 

8^9 (1579) 

871 (1599) 

Working fluid inlet 
temperature, *C (*F) 

580 (1076) 

580 (1076) 

579 (1075) 

579 (1075) 

Working fluid out let 
temperature, *C (°F) 

756 (1393) 

756 (1392) 

754 (1390) 

754 (1390) 

Maximum HSA outer 
shel 1 temperature, 
*C CF) 

81 (178) 

81 (178) 

83 (181) 

8? (182) 

AP/P, percent 

0.504 

0.504 

0.504 

0.504 

HSA heat loss, percent 

5.84 

5.90 

6.17 

6.29 


In order to increase the PICS temperatures to comply with the specified 
requirements, the thermal resistance between the PICS and the working fluid must 
be increased. One way to do this is to decrease the surface emissivity of the 
Hastel loy X HSHX. This would seem to be a reasonable approach, since to attain 
the baseline design emissivity of 0.8, surface treatment is required. The 
stability of this high emissivity is problematical. 
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A sonsitlvity study was conducted by performing HSA thermal analyses for 
HSHX's of various surface emlsslvitles. It should be noted that uniform emls- 
s. vities w«-e assumed for the entire outside surface of the HSHX — the portion 
viewing the IHS and the portion viewing the HSIS. the emlsslvlty of the surface 
of the HSHX facing the HSIS does not directly affect the thermal resistance 
between the PICS and the working fluid, but does affect the HSHX temperature. 

The results of the study are summarized in the remainder of table 17.^ Note 
that for each case, the maximum HSHX wall temperature was 788*C <1450*F). 

It can be seen that the 0.25 emisslvity yields PICS temperatures very 
close to the required temperature range of 1150* to 1200*C. This emissivity 
corresponds to the value for untreated Hastel loy X. Naturally, the temperature 
level of the entire IHS increases as the PICS temperature increases, as can be 
seen from the ablation sleeve temperatures. The decreased emissivity does cause 
a small performance penalty, however. The 0.25 emissivity case results in a 2 C 
(3*F) drop in the attainable working fluid outlet temperature. In addition, the 
heat loss from the HSA is increased tor the lower HSHX emissivity. 

Based on the results of the sensitivity study, the untreated Kfc>HX 
enissivity of 0.25 looks promising; however, the effect of this emissivity must 
be checked for all operating conditions Including the ECS, ACS, and the startup 
transi ents. 

As previously mentioned, the HSA heat sink assumed for the preceding 
analyses was 49*C (120*F). To estimate the effect of the heat sink on the^ 
Important HSA parameters, an analysis was performed for a heat sink of -62°C 
(-80*F ) . The change In heat sink temperature had a minimal effect on the 
key HSA parameters; the temperatures of the PICS and HSHX decreased by less 
than 0.3*C (0.5*F) and the HSA heat loss Increased by about 1 W. This would 
seem to Indicate that diurnal temperature changes in orbit should have a 
negligible effect on HSHX thermal fatigue. 

Auxiliary Cooling System 

The ACS is required to provide cooling for the HSA during all nonoperat iona I 
phases of the Brayton cycle when there is no circulating working fluid available 
to remove The heat produced by the IHS. Thus, the ACS must operate on the 
launch pad, in space prior to system startup, and during any subsequent system 
shutdown. The requirements for ACS operation are as follows: 

(1) On the launch pad, all exposed surfaces must be maintained below 
193*C (380*F) to prevent possible Ignition of fuel vapors present 
In the prelaunch environment. 

(2) The PICS temperatures should be maintained in the 1050* to 120C*C 
range — the launch and ascent requirements given In table 1. It 
should be noted that the PICS temperature "window" is wider for 
ACS operation than for normal operating conditions. 

(3) The HSHX temperature should not be so high as to seriously degrade 
fatigue life. 
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dome to below 193*C (380 # F). The cocoon Is much cooler then the end dome and 
does not require additional cooling. 

In order to determine the appropriate fill gas, thermal analyses were 
conducted using helium, neon, argon, and xenon. The baseline design HSHX, with 
a surface emlssivity of 0.8, was assumed for the fill gas comparison study. A 
sunnary of the important HSA temperatures for the different fill gases is given 
In table 18. The xenon results are not shown since with this fill gas the 
HSA temperatures are so high that melting of some of the Insu lotion foi 1: s was 
predicted. With helium as the fill gas, additional cooling of the HSA top end 
dome is not required. 

TABLE 18 

HSA TEMPERATURES FOR LAUNCH PAD AUXILIARY COOLING OPERATION 
I ; Temperature, °C (°F) 




s I eeve 
(maximum) 

HSHX 

(maximum) 

Top end 
dome 

(maximum) 

Cocoon 


851 

(1564 

to 897 
to 1647) 

1014 

(1858 

to 1050 
to 1922) 

1219 to 1249 
(2227 to 2280) 

1051 to 1085 
(1923 to 1985) 

619 

(1147) 

782 

( 1440) 

1023 (1874) 

837 (1538) 

573 

(1064) 

751 

(1383) 

1007 (1844) 

771 (1419) 

183 

(362) 

152 

(305) 

161 

(322) 

153 

(307) 

77 

(171) 

74 

(166) 

74 

(165) 

74 

(166) 


a No additional cooling for top end dome. 

^HSHX surface emissivlty of 0.25. 

It can be seen from the first three columns of table 18 thal^the PICS 
temperatures come closest to the specified range of 1050° to 1200°C with 
neon as the fill gas. Helium provides surh a good heat transfer path that the 
predicted PICS temperatures are far below the desired range. With argon, the 
heat transfer resistance is not sufficiently decreased to cool the P CS below 
1200*C. In addition, the HSHX maximum wall temperature using ar 9°n # ' s ~ 
<394°F) above the maximum steady-state operating temperature of 788 C (M50 M 
The variation In PICS tanperature Is in the HSA axial direction, and is due to 
heat transfer from the ends of the HSA. The absence of fluid flow In the HSHX 
allows for more uniform PICS temperatures than for normal operating conditions 







The complete HSA temperature map 
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TABLE 19 


HSA TEMPERATURES FOR SPACE AUXILIARY COOLING OPERATION 
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pU fl + , CO ' T ’ pa^e<, " ith +he 49 *C (120-F) heat sink, the maximum 

PICS temperature increased by about 3*C and the maximum HSHX temperature 
increased by about 4*C (7*F). This Indicates that the heat sink temperature 

°I. ACS opera+,on tha " on normal steady-state operation. 

J CS ? pera | ,on ’ the sensitivity of the key HSA temperatures to the 
heat sink temperature Is not large enough to have a significant system impact. 
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The fueling Is Initiated by removing the EHS from the BIPS HSA and 
replacing it with the IHS. Since this operation Is likely to take a reasonable 
amount of time, the IHS can be assumed to reach a steady state with respect to 
the neon environment. These steady-state IHS temperatures are the Initial 
temperatures for the fueling transient. 

To determine the steady-state IHS temperatures, a thermal analysis was 
conducted using the HSA thermal model depicted in fig. 15. Only the IHS nodes 
were used. The heat sink for the IHS was assumed to be stagnant neon at ambient 
pressure and 27°C (80°F). The results of this analysis are summarized in table 
20 . 


TABLE 20 

AVERAGED IHS STEADY-STATE TEMPERATURES IN THE MLS 


Component 

Temperature 

°C 

•F 

PICS 

791 

1455 

GIS Inner surface 

721 

1329 

GIS outer surface 

662 

1223 

End crush-ups 

495 

923 

Retaining ring 

397 

747 

Aeroshel 1 

337 

638 

Aeroshel 1 end caps 

326 

618 

Ablation sleeve 

279 

535 


The IHS temperatures shown in table 20 are well below the temperatures 
at normal ACS or BIPS operation. This is due to the greatly reduced thermal 
resistance when the IHS is outside of the HSA. Note that the temperature of 
the outer surface of the IHS is only 279°C (535*F). 

Since the fueling transient takes a few hours to approach steady state, a 
simplified HSA model was required to reduce computer time. In the simplified 
model, each set of axial components in the IHS, HSHX, HSIS, and HSA outer case 
Is represented as a single node. In addition, the HSHX quadrant nodes in the 
coil circumferential direction are lumped together as a single node. 

The simplified model has one PICS node and four HSHX nodes. Since there 
is no fluid flow through the HSHX during the fueling transient, the preceding 
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assumptions are not unreasonable. However, due to heat losses out the ends of 
the HSA, there will be a small axial temperature variation. The PICS and HSHX 
temperatures predicted with the simplified model represent the average axial 
temperatures. The actual maximum PICS and HSHX temperatures will thus be some- 
what higher than the average temperatures. 

The temperatures listed in fable 20 are input as the initial IHS tempera- 
tures. All other initial HSA temperatures (including the HSHX) as well as the 
constant heat sink temperature are assumed to be 49°C ( 1 20 * F ) • The HSHX is 
assumed to be filled with stagnant Xe-He working fluid at 0.21 MPa (30 psia), 
the approximate BIPS pressure for zero-flow conditions. 

The predicted HSHX temperatures during the fueling transient are shown in 
fig. 20 . Also shown are the PICS temperatures. The reference HSHX, with an 
emissivity of 0 . 8 , was used in the model. It can be seen that after the initial 
rapid HSHX temperature rise, the approach to steady state is slow. The circum- 
ferential temperature gradients in the HSHX tube do not seem to be particularly 
severe during the transient. The predicted PICS temperature is around 1030°C 
after 6 hours. 


BIPS Startup Transient 

When cool working fluid enters the HSHX during the system startup, large 
thermal gradients may be produced in the heat exchanger, producing stresses 
that must be evaluated to determine the LCF life of the unit. The HSA has a 
requirement for a minimum of 100 startup cycles without evidence of fatigue 
cracking. The study work concentrated on the thermal gradients and cycle life 
prediction for the HSHX. The thermal analysis is discussed in this section. 

The LCF life predict’on is discussed in tho structural analysis section. 

The BIPS program office at AIRRHX has provided the latest BIPS startup con- 
ditions for the proposed flight system start (B. Minshall, BIPS-GDS Coordination 
Memo A0609, April 21, 1978). The HSHX inlet flow and temperature from rhe AIRPHX 
system transient computer program are shown for the BIPS startup in fig. 21. The 
conditions are identified as Run 600. Also shown are the HSHX outlet tempera- 
tures as predicted from the system program and the detailed HSA model. These 
temperatures will be discussed later. Note that the flow is for the entire 
two-HSA BIPS; the flow through each HSA is taken as one-half of the total. It 
should be emphasized that the startup conditions are preliminary as the flight 
system scenario has not been finalized. 

The proposed startup scenario begins with the HSA In the steady-stale 
space ACS mode of operation. The HSA heat sink is assumed to be about 16°C 
(60 *F). This is expected to be the coldest heat sink (radiator temperature) 
for the HSA while the BIPS is in the spacecraft; the coldest heat sir.k repre- 
sents the most critical startup condition. The BIPS rotating equipment will 
be motored and allowed to accelerate to full speed to check engine operation. 

At this point, the neon is vented from the HSA. This is the zero time for the 
transients shown in fig. 21. it can be seen from fig. 19 that the stagnant 
working fluid in the HSHX during the ACS deployment is close to normal operating 
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temperatures. The scenario Is described In detail In the BIPS memo referenced 
In the preceding paragraph. 

To analyze the startup temperatures, the thermal analyzer computer program 
was run in the transient mode using the complete HSA model shown in fig. 15. 

The initial conditions for the transient were taken from a space ACS analysis, 
similar to that shown in fig. 19, but with a heat sink temperature of 16°C 
(60 # F) . The temperatures from this run are very similar to those shown in 
fig. 19. 


The key results of the transient analysis are shown in fig. 22. HSHX 
temperatures are shown at the inlet and outlet ends of the unit. The reference 
HSHX, with an emissivity of 0.8, was used In the model. Temperatures for two 
of the four circumferential nodes are shown; the node facing the IHS and the 
node facing the HSIS. Complete HSA temperature profiles are available from 
the author for various times during the startup transient; these are not 
presented in this document. 

It can be seen that the Inlet portion of the HSHX, Initially hot during 
ACS operation, Is rapidly cooled by the cold working fluid. Since tho IHS 
remains hot, the heat flux to the inlet region of the HSHX increases, resultina 
in an increase in tube circumferential temperature gradient. The inlet end 
minimum tube wall temperature and maximum circumferential temperature gradient 
occur about 180 sec after startup. The temperatures at the HSHX outlet end are 
still decreasing at 420 sec; these temperatures must bottom out and increase as 
steady state is approached. 

Referring back to fig. 21, the predicted fluid outlet temperature is 
/shown and compared with the HSHX outlet temperature from Run 600. The AIRPHX 
system computer program has a simplified model for the HSA and HSHX; thus the 
outlet temperatures will be somewhat different than for the detailed HSA model. 
The agreement must be close enough to ensure that the HSHX inlet temperatures 
are reasonably accurate. The agreement in fig. 21 would appear to be satis- 
factory for this analysis; for a more accurate thermal analysis, another itera- 
tion between the system, model and the HSA model would be required. 

The PICS temperatures are Initially at the space ACS values. During the 
startup, the PICS temperatures increase and approach the steady-state values. 

The steady-state temperatures are not exceeded during the startup transient. 

Emergency Cooling System 

The iCS is required to maintain the PICS at a safe temperature (see table 1) 
In the evont of a failure In the BIPS, such as a loss of working fluid accident 
or a breakdown in the ACS. The ECS must remove the heat generated by the IHS 
fuel that is normally removed by the circulating working fluid. 

The present ECS design relies on the successive melting of the HSIS foil 
layers. If a system failure occurs, the heat released by the fuel is Initially 
absorbed by the HSA components, including the IHS, the HSHX, and the HSIS. 

Only a small amount of heat can be rejected to the environment due to the high 
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:ure 22.-- HS‘ Temperature During BIPS Startup. 



thermal resistance of the Insulation. The absorbed heat causes the temperature 
of the HSA components to rise. As the temperature of the hottest foil (closest 
to the IHS ) reaches the foil melting temperature, the layer melts and allows 
increased irradiation of the next foil, with a subsequent increase in tempera- 
ture. As each foil melts, the resistive path through the Insulation bundle 
decreases, and the heat loss to the heat sink Increases. At some point, a 
steady state will be reached, with a certain number of foil layers surviving. 

At this point, the heat generated by the fuel equals the heat rejected to the 
environment. As listed in table 1, the PICS temperature must be controlled 
both during the transient meltdown and the post-meltdown steady state. 

With this goal in mind, both transient and steady-state thermal analyses 
were performed tor the ECS mode of operation. Since the ECS mode represents 
a long running transient, the simplified HSA model was used. In this model, 
as discussed previously, each set of axial components is represented as 
a single node. Use of the simplified model results in predicted average PICS 
temperatures somewhat below the maximum PICS temperatures. 

In order to determine the transient meltdown characteristics, it was 
necessary to modify the AIHLA thermal analyzer program to account for foil 
melting. To reduce computation time, the 60 foil layers in the HSIS were 
lumped Into 14 nodes with the innermost node composed of 3 foils, the next 
eleven nodes five foils each, and the last two nodes one foil each. The 
foil nodes or packets are modeled to melt when the node reaches the melting 
temperature. The melted packet is assumed to resolidify on the next outer 
packet. The mass of the melted packet is combined with the outer packet to 
form a new node of increased thermal capacitance. This reduces the effective 
number of radiation shields and increases the thermal conductance through the 
HSIS. At any time, the innermost node contains. In addition to its own mass, 
the mass of all the previously melted foils. In this manner, the heat of 
fusion required to melt the foils is recovered upon resolidification. 

For the present analysis, all of the foils were assumed to be a mixture 
of 40 percent Ni and 60 percent Pd, with a constant melting temperature of 
1237*C (2259*F). There seemed to be no reason to include any higher melting 
point pure Ni or Mo foils, since it became apparent that it would be desirable 
to melt as many of the foil layers as possible. 

The TECO correlation tor heat transfer through the foils, eg. (1), was 
not used tor the ECS transient since it represents heat transfer from the hot 
side to the cold side of the entire toil package. Due to the change in effec- 
tive foil emlsslvlty with temperature, eq. (1) is not valid tor each of the 
insulation nodes used for the ECS iranslent analysis. Instead, eftecrive foil 
emlsslvlty data, presented in T. Ashe, 8IPS-G0S Coordination Memo A0090, December 
23, 1975, were used so that each node would have an emlssivlty appropriate to its 
temperature. The effective emissivity includes both the radiation and conduction 
heat transfer in the toils, modeled as radiation heat transfer; the conductive 
path Is provided by the refractory oxide binder in the Multi-Foil system. The 
emissivity data were used with the following equation for radiation heat transfer 
■through a series of radiation shields (ref. 7). 
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K-.*] 


whore 4> Is the heat flux perpendicular to the foils 
o is the Stefan-Boltzmann constant 

e f is the effective foil emisslvity at the average temperature of a 
packet 

N is the number of foil layers In the packet 
T h is the hot-side temperature of the packet 
T c is the cold- side temperature of the packet 

During the ECS transient, the Hastel loy X HSHX and insulation support 
mandrel are heated well above their normal operating temperature.. Indeed, 
during the transient the temperatures exceed the solidus point of Hastel loy X, 
1260*C (2300 °F ) . in the meltdown analysis, it was assumed that the Hasteiioy* 

X components would r<«nain in the model as thermal barriers unless their tem- 
peratures exceeded :h.c- liquidus point, 1354°C (2470°F) 

Other assumptions made for the ECS transient analysis Included the 
fol lowing: 

CD The heat sink was 49*C M20*F). 

(2) The reference design HSHX was used (surface emisslvity = 0.8). 

(3) Any effects due to sublimation of either the insulation or the HSHX 
were not included in the model. 

(4) The starting point for the ECS transient was the steady-state normal 
operating condition. 

The key results of the ECS meltdown transient are shown in fig. 23. It 
can be seen that during the transient, the HSHX t«nperature does not exceed 
the liquidus point of Hastel loy X. Thus, the HSHX and the support mandrel 
(which is cooler than the HSHX) are retained In the model as thermal barriers 
throughout the meltdown. Elimination of these nodes would tend to decrease 
the HSA resistance and ultimate PICS tanperature. 

Starting at Its steady-state temperature at zero time, the PICS tempera- 
ture Increases until the first foil packet is melted at about 130 min. At this 
point the PICS temperature is 1432*C. Foil meltdown continues for about 40 min. 
At this point (as can be seen from the heat loss curve in fig. 23) the thermal 
resistance of the HSIS has decreased sufficiently so that almost all of the 
heat produced by the foil is rejected to the environment. The remaining 
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Figure 23. — ECS Meltdown Results. 
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foils cannot continue to increase in temperature so the meltdown Is terminated. 
The analysis predicts that 12 cylindrical and 7 planar (end) foils remain. Due 
to the lumping of foil layers Into packets, the next unmelted nodes contain 5 
toil layers each. Thus one more melted cylindrical node would result In 7 layers 
remaining and one more melted planar node would result in 2 layers remaining. 

The PICS temperature reaches a maximum of 1474*C at about 175 min. Since 
the HSA is close to steady state, the PICS temperature Is only 16 # C above its 
ultimate steady-state temperature of 1458*C. This is considerably higher than 
the PICS temperature limitation of 1200*C for ECS post-meltdown steady state. 

The absolute temperature limit during the transient of I600 # C is not approached 
during the meltdown. 

A transient analysis was performed replacing the inner Hastel loy X insula- 
tion mandrel with a mandrel made of a material similar to the main foils (40J 
Ni-<>0f Pd). This sheet would melt during ECS operation and reduce the thermal 
resistance. The results Indicated that 12 cylindrical and 7 planar foils 
remained after meltdown — exactly the same as with the Hastel loy X mandrel. Due 
to the removal of the mandrel as a radiation shield, the steady-state post- 
meltdown PICS temperature was lower by about 10*C. 

To assess the influence of the heat sink temperature on the ECS operation, 
several analyses were performed using a sink temperature of -62*C (-80°F). 

Very little change was noticed for either the transient or the ultimate steady- 
state condi tions. 

Post-mel tdown steady-state analyses were performed for various numbers 
of remaining foil layers, without regard to how the sfeady stato was achi .e'. 

To conserve computation time, the simplified HSA model was used. The PICS 
and insulation temperatures for steady state are shown in table 21. It can 
be seen that to maintain the PICS temperature below the specified limit of 
1200°C (table 1) all of the foils must pe removed. 


TABLE 21 

ECS POST-MELTDOWN STEADY STATE 


Remain 1 ng foils 

PiCS 
temp. , 

•c 

Inner toil temp., °C 

Cy 1 1 ndr ica 1 

Planar 

0 

1160 

- 

- 

1 

121 1 

645 

b?2 

2 

1286 

878 

804 
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The effect of foil sublimation on the ECS operation has not been analyzed. 

rra r -. 

dr-oD Table 21 shows the temperatures of the innermost foil with one and 
Ten iers remaining. These temperatures are .ell below the melt.ng 
point (1237°C) and the vapor pressure would be quite low. t seems very 
unlikely that all of the foils remaining after meltdown would be removed by 
sublimation, especially when giving consideration to the 

of the melted foils. The sublimation phenomenon is certain y an area 
that requires further investigation given the necessity of removing a of 
the toils to meet the PICS temperature requirements. The possib e sublima 
tion of the Hastel loy X HSHX and/or insulation support mandrel also shou d 

be investigated. 

To attain the desired temperature level with the present type of ECS 
(i.e., foil meltdown) would probably require using HSIS foil material wit 
a lower melting point than the Ni-Pd. With a lower melting point it is 
likely there would be a higher vapor pressure; this must not be high enough 
to endanger foil stability under normal operating conditions. It may a so 
be helpful to reduce the number of foil layers used. This, o course would 
tend to increase the heat leak under normal operation. The effect of such 
^change on the ACS would also have to be considered. The use of a fusible 
inner insulation mandrel should be considered as a design option. 

The ECS analysis just discussed was for the reference design HSHX with 
an emissivity of 0.8. It appears that an untreated emissivity of 0.25 may 
be more desirable for the other modes of operation. This lower emissivity 
should be included in future ECS studies. Since the reduced emissivity 
increases the HSA thermal resistance, the PICS temperature problem sou 
become more difficult unless the increase in HSA temperature causes 
HSHX to exceed the Hastel loy X liquidus point. If this occurs, the HSH 
is relieved as a thermal barrier, reducing the HSA thermal resistance. 

Results of the recent TECO -mini -meltdown" tests for Ni-Pd Multi-Foil 
insulation have become available (T. Ashe, BIPS-GDS Coord . ant ion Memo A0596 
March 6, 1978). These results are essentially consistent with the transient 
meltdown analysis discusseo in this section. 
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STRUCTURAL ANALYSIS DESCRIPTION AND DISCUSSION OF RESULTS 


The structural analysis performed during the study concentrated on satisfy- 
ing the following requirements: 

Pressure Containment— The HSHX Is designed for pressure containment 
for a 7-yr life. A maximum of 1 percent creep deformation Is allowed 
for the required life. 

Low-Cycle Fatigue — The HSA must withstand a minimum of 100 startup 
cycles. The analysis was limited to the HSHX. 

R esonant Frequency — In accordance with the BIPS requirement, no HSA 
natural frequencies should be below 60 Hz to avoid any spacecraft- 
induced resonances. 

Dynamic Loadinq— The HSA must withstand launch loads of shock, vibra- 
tion, and acceleration without any loss in structural integrity. 

M IL-STD-8 10B is the specified requirement for dynamic loading. An 
equivalent static loading of +50 g, where g is the normal acceleration 
of qravity, was used In the analysis and compared with MIL— STD-810B. 

In the calculation of structural properties for the HSHX, a design margin 
ot 28 °C (50°F) was used. Thus, for a HSHX temperature of 788°C (M50*F), struc- 
tural properties at 816°C (1500°r) are used. 

Pressure Containment 

The HSHX is designed tor pressure containment tor a 7-yr life. The design 
allowable stress is based on 1 percent creep deformation, at the appropriate 
fluid pressure level, accumulated over the life of the unit. 

For short-term stress, the design practice employed by AiResearch is to 
size the pressure-carrying structure for proof pressures of 1.5 times the 
working pressures and for burst pressures of 2.5 times the working pressures. 

The structure must not yield at proof pressure or rupture at burst pressure. 

This Implies that the proof pressure is the governing design condition if the 
ratio of yield stress to ultimate stress Is less than 0.6 and that the burst 
pressure will govern if the ratio Is greater than 0.6. 

When the limiting stress Is due to bending, a small amount of yielding 
can be allowed in the outermost fibers, which leads to a modified stress distri- 
bution through the thickness. Tho Ideal plastic bending moment is 1.5 times 
the computed elastic bending moment for the same peak stress. Accordingly, 
the allowable Indicated elastic stress due to bending loads is taken to be 
1.5 times the allowable values for uniform tension or compression loads. 
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At the elevated operating temperatures and long iife requirement for the 
HSHX, if has been found that the creep stress is always the controlling stress. 
The 1 percent creep deformation criterion was used to determine the required 
wall thicknesses for all pressure containment components in the HSHX — tubes, 
elbows, and manifolds. An aging factor, F a , and a tube weld efficiency, n. 
were used to reduce the allowable stress for all pressure containment designs 
as fol lows: 


CT a 



(14) 


where o a is the allowable stress and 0 C is the 1 percent creep stress for 7 yr. 

The 1 percent creep stress is taken at the maximum temperature of the component 
including the 28 °C (50°F) design margin. Potential degradation of material 
properties due to contamination, radiation- induced swelling, and radiation- 
inouced creep have not been considered. An aging factor of 1.2 was assumed 
for a 1 1 designs. 

Configuration study . — For the configuration studies the long-term creep 
stresses were estimated from refs. 3 and 9 for the C-103 and the Haste I Icy X 
designs, respectively. A tube weld efficiency of 0.95 was assumed for the 
tubes and manifolds, a resonable value for rolled and planished butt welds. 

I be required tube and manifold wall thickness, t, is calculated by the thin- 
wal I pressure containment relationship: 



(15) 


where K is a tube wall thinning factor 

P is the design pressure containment criterion, 0.793 MPa (115 psia) 

D is the tube OD 

For the single-tube helical coil, the helical-flow tube bank, and the 
toroidal manifolds for both multitube designs, a wall thinning factor of 1.05 
was assumed to account for any thinning due to tube bending. The dimpled 
tubes of the axial-flow tube bank design are not bent, but the dimple forma- 
tion results in some wall thinning in Ihe dimple area. For the relatively 
shallow dimples used in the design, AiResearch experience indicates a maximum 
wall thinning of around 12 percent. Thus, a K of 1.12 was used for the dimpled 
tube wall thickness determination. 

Selected Configuration . — Required tube and elbow wall thicknesses were 
determined for the selected single-tube helical-coil Hastel loy X configuration 
for coils of various numbers of turns (see table 12). Updated creep stress pro- 
perties for Hastel loy X were used (J. Hadley, BIPS-GDS Coordination Memo A0521, 
September 21, 1977). The 1 percent creep stress at 816°C ( 1 500°F )— 788°C 
( 1 4 50 0 F ) maximum tube wall temperature plus the 28°C (50°F) design margin— is 
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15.2 MPa (2200 psi). In the determination of the allowable stress In eq. (14), 
a tube weld efficiency factor of 0.95 was assumed for the HSHX tube. In the 
transition elbows, which connect the tube with the BIPS ducts, the radius of 
curvative precludes roll planishing, so a lower weld efficiency of 0.7 was 
assumed. The Inlet and outlet elbows were designed to have identical wall 
thicknesses even though the inlet elbow is at a considerably lower temperature 
(and hence, higher allowable stress level) than the outlet elbow. The tube wall 
thicknesses were calculated using eq. (15) for the design pressure containment 
criterion of 0.703 MPa (102 psia) and a tube wall thinning factor due to bending 
of 1.05. 

The results of these calculations are shown in table 22. As previously 
discussed, an additional allowance of 5 percent was added to the minimum 
required tube wall thickness to account for manufacturing tolerances. This 
allowance is reflected In the nominal wall thicknesses listed in table 12. 
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TABLE 22 

TUBE WALL THICKNESS 


i 


H*J 


i: 


SB? 
- ,1 


Number 
of turns 
in coi 1 

Tube 00, 
cm (in.) 

Mini mum 
tube wal 1 
thickness, 
mm ( In. ) 

Minimum 
elbow wal 1 
th ickness, 
mm (in.) 

9 

4.42 (1.74) 

1.35 (0.053) 

2.36 (0.093) 

10 

3.96 (1.56) 

1.22 (0.048) 

2.11 (0.083) 

11 

3.61 (1.42) 

1.12 (0.044) 

1.93 (0.076) 

12 

3.30 (1.30) 

1.02 (0.040) 

1.75 (0.069) 

13 

3.05 (1.20) 

0.94 (0.037) 

1.63 (0.064) 


As a check on the selected 10-turn design, a three-dimensional thin-shell 
finite-element model was constructed and the pressure containment stress 
calculated using the ANSYS computer program. The ANSYS is a general-purpose 
finite-element computer program used for stress and vibration analysis. The 
program is discussed in appendix C. The model was for a single turn of the 
HSHX with boundary conditions applied to simulate the effects of the continuous 
tube. A model also was constructed for the transition elbows. The tube and 
elbow models are shown in figs. 24 and 25, respectively. 

The results of the ANSYS analysis substantiated the 10-turn minimum wall 
thickness calculation and verified the applicability of eq. (15). 
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Thermal Stresses and Low-Cycle Fatigue 

The HSA has a requirement for a minimum of 100 startup cycles without 
evidence of fatigue cracking. The present analysis concentrated on the ICF 
life prediction of the HSHX. The temperature profiles calculated from BIPS 
startup transient Run 600 were used as inputs to the thermal stress computa- 
tions (see the thermal analysis section). 

A finite-element model of the HSHX was constructed for use with the ANSYS 
computer program to compute the thermal stresses. The HSHX has a continuous 
strip that slips into the separation between successive turns of the coil and 
is brazed to the coil (see the system design section). The brazed strip is 
implemented in the ANSYS model by the use of appropriate boundary conditions 
at the required angular location of each turn in the coil. The brazed joint 
Is assumed to have a braze efficiency of 0.5. The model Is shown in fig. 26. 

Thermal stresses are set up in both the axial (flow) and circumferential 
directions around the tube. The steady-state thermal stresses are relatively 
low, well below the yield stress. With the expected stress relaxation, it 
appears that the steady-state thermal stresses will not impact on the required 
7-yr life. 

The transient analysis indicates that the most critically stressed region 
of the HSHX is the inlet to the tube. This is due primarily to the large 
circumferential temperature gradients resulting from the cold inlet fluid 
temperature (see the thermal analysis section). It can be seen from fig. 22 
that the largest temperature gradients at the inlet end occur around 180 sec. 
Although the yield stress is lower at the outlet (hot) end, the smaller tem- 
erarure gradients result in less critically stressed regions. The axial 
temperature gradients produce much smaller stresses; the maximum axial thermal 
stress is only about one-seventh of the maximum circumferential thermal stress. 

The maximum combined thermal stress due to circumferential and axial 
temperature gradients is only about 50 percent of the yield stress at the Inlet 
end. In the absence of stress concentration effects, there would be no plastic 
strains generated. However, the LCF life of the HSHX will be I imited by the 
repetitive creep relaxation; to account for the relaxation a stress concentra- 
tion factor of 3 was assumed. This large value was assumed so as to include 
the effects of the local changes in coil thickness produced hy the continuous 
braze strip and the HSHX mount brackets. 

Based on a comparison with a fatigue life study conducted under a DOE 
contract (ref. 10), an l£F life of 570 startups Is predicted. This is 
well above the HSA requirement of 100 cycles. The faTlgue failure location 
Is at the braze Joint at the cold end of the tube. It should be noted, 
however, that only the startup portion (Run 600) of a complete cycle was 
considered. Reductions in fatigue life due to other modes of operation 
(e.g., steady state, HSA fueling transient, etc.) have not been included. 

These effects are expected to be small. 
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Resonant Frequency 

The BIPS has the requirement to maintain ^”t?™^uSeUs!“Tl.ls 

above 60 Hz to avoid '7 t P °U' “HSHX° U «al the main driving force in the selection 
requirement, as applied to the HSHX, , , as ttle natural frequency of 

of the continuous brazed s rip T he brazed strip has the additional advantages of 

the free coil Is below 60 Hz. TheM * on effect of me internal pressure, 

taking up the un.lndlng load tueto the tube, and helping to mainta n 

significantly reducing tetec^o of me~nt I nuous Prized strip coll design is 
direct alignment. Ti e ° t his study concentrated on determining 

discussed in ref. 13. Tn. Amounts. ensuring that all modes 

£eTt r he ^STeSS-^^'ng these frequencies -Itb the calculated 
natural frequencies of the HSHX. 

. i t uquy was constructed to determine 

A finite-el^ compute, to that shown in 
the dynamic response of the un. that attach the heat exchanger coil 

fig. 26. The forwa H r ^ nd a / f / hr ^! C dimensio nal beams with stiffnesses equ. valent 
to the HSA were modeled as ihr * assumed to be held in the radial 

to the brackets. The beam elemen . a+ the aft graphite plate. The 

and axial directions but were re ^ freely in the' fore and aft direction 

forward beam elements were a °" |n +he radia i directions by the attachment 
and to rotate freely, but we The coj t itself was modeled using 

to the forward graphite plate In the • ^ co J|-to-coil brazed strips were 

three-dimensional curved tubular e^em© 0 ^ elefTie nts. The axial and bending 

set to yield tbe .eight of the brazed strip. 

_ ANSYS computer program to determine the HSHX 

resort me th r e es A u| S ts S ,or ?be first fe. nature, frequencies ere 

gi ven in tab I e 23. 


TABLE 23 

HSHX NATURAL FREQUENCIES 


Frequency, 

Hz 


Lateral bending mode (bending of coil) 

Axial mode (resonance of mount brackets) 

Lateral rocking mode (side displacement of brackets) 
Undrf ined higher order modes — — 
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The first and third nodes are Illustrated in figs. 27 and 28, respec- 
tively. As can be seen from table 23, all natural frequencies of the HSHX 
are wel I above 60 Hz. 

A finite-element model of the entire HSA was constructed to determine the 
HSA natural frequencies at the HSHX mount points. The model representation 
and restraints are sunmarlzed In table 24. 


TABLE 24 

HSA DYNAMIC MODEL REPRESENTATION AND RESTRAINTS 


HSA feature 

Model representation 

Model restraint 

Forward and aft 
end domes 

Shel 1 el emerts 

Center of aft dome fixed In X-Y 
(lateral) plane to simulate mount. 

Outer shel 1 

Shel 1 el ements 

Outer shell, forward dome, and 
spider fixed in X-Y-Z directions 
at three locations to simulate 
ball joint suspension. 

Forward and aft 
spiders (end 
enc losures) 

Centra 1 hub as shel 1 
element joined to 
outer shel 1 with six 
3-D beam el ements 

Beams coupled to outer shell in 
X-Y-Z directions but free to 
rotate. 

Preload screws 

Shel 1 el ements in 
conical shape 

Rigid attachment to spider; 12 
compression-only spar elements at 
graphite shoe Interface (Z direc- 
tion). Coupled to graphite shoe 
in X-Y plane (lateral restraint). 

Graphite shoes 

24 shell elements 

Graphite shoe - IHS interface: 12 

compression-only spar elements 
restrained in X-Y plane to repre- 
sent socket design. 

IHS 

12 shell elements 
mod i f i ed to represent 
actual weight and 
sti ffness 


HSHX 

8 shel 1 el ements 
modified to represent 
actual weight and 
sti ffness 

Aft graphite shoe: coupled in the 

X-Y-Z directions. 

Forward graphite shoe: coupled 

in the X-Y plane only. 
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The model was run on the ANSYS computer program to determine the HSA 
natural frequencies at the attachment points of the HSHX to the aft (outlet) 
graphite support plates. The frequencies are listed in table 25. 


TABLE 25 

HSA NATURAL FREQUENCIES AT FT END 


Mode 

Frequency, 

Hz 

Torsional mode (central mass rotating with legs of 

28 

spiders acting as springs) 


Axial mode (central mass moving axially with legs of 

92 

spider acting as springs) 


Lateral rc :ing mode (central mass moving laterally with 

330 

trunnion action of preload screws causing spide r to 


di stort) 


Undefined higher order modes 

754 


The lateral mode is illustrated in fig. 29. 

Th*' torsional mode natural frequency is below the design goal of 60 Hz. 
In addit.on, a comparison of -fables 23 and 25 indicates the possibility of 
mode excitation of the axial mode of the HSHX (342 Hz) by the lateral rocking 
mode of tne HSA (330 Hz). This could occur even though the two mode shapes 
are different. The fourth modes of the HSA (754 Hz) and the HSHX (734 Hz) 
are close together; however, these modes are complex and will be difficult to 
excite. In addition, these frequencies are above the anticipated sinusoidal 
excitation range for the HSA of 2 to 500 Hz. No coupling seems possible 
between any other modes. 

To increase the natural frequency of the torsional mode of the HSA, it 
is proposed to add material to the spokes of the spiders (end enclosures) to 
increase their lateral stiffness. Analysis indicated that increasing the 
average width of the spokes from 1.63 cm (0.64 in.) to 2.64 cm (1.04 in.) 
would be sufficient to raise the natural frequency to 60 Hz. The additional 
material also would increase the stiffness of the spokes, and hence natural 
frequency, in the axial and lateral directions. The HSA natural frequencies 
for the modified end enclosures are presented in table 26. The unmodified 
HSA natural frequencies from table 25 and the HSHX natural frequencies from 
table 23 are included for comparison. 
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Figure 29.— HSA Lateral Rocking Mode. 




J I L__I J 1 I L— i I 1 L I I I — • 

TABLE 26 

NATURAL FREQUENCY SUMMARY 


HSA 

HSHX 

Frequency, Hz 


Frequency 

Hz, 


Original 

Mod if 

ed 

Vibration mode 

Vibration mode 

28 

60 


Torsiona 1 

261 

Lateral (bending) 

92 

118 


Axial 

342 

Axial 

330 

420 


Lateral (rocking) 

648 

Lateral (rocking) 


Since the resonant frequency of the lateral rocking vibration mode of 
the HSA increased to 420 Hz, the possibility of coupling between the HSA 
and the HSHX has been eliminated. It is recommended that the end enclosure 
modification be made to the HSA design. 

Dynamic Loading 


The HSA must withstand launch loads of shock, vibration, and acceleration 
without any loss in structural integrity. The requirement for dynamic loading 
is MIL-STD-810B, except that the damping ratio is specified as 5 percent. To 
simplify the computations, an equivalent static loading in terms of g levels 
was assumed. A level of j+50 g was used. This level was shown to be well above 
the maximum response from the MII-STD-810B inputs. The temperatures used for 
the determination of the allowable stresses were those just prior to launch, 
i.e., the ACS temperatures shown in fig. 18. The 28°C (50°F) design margin 
was app I i ed. 

The finite-element model of the HSHX (see fig. 26) was used to determine 
the loads at the HSHX mount points due to the application of 50-g static loads 
in the axial (+Z) and two radial directions (+X and +Y). The aft (ho+ end) mount 
points for the HSHX were assumed to be restrained in the axial, radial, and 
tangential directions. The forward (cold end) mount points were assumed to be 
restrained in the radial and tangential directions only. The loads at the aft 
end were found to be greater than those at the forward end. The calculated 
loads are presented in table 27. 






These loads were app I i ec to the HSHX mount brackets and to the combined 
IHS and HSHX support plate in the HSA finite-element model (see table 24) to 
determine the structural integrity for the anticipated dynamic loading. For 
the support plate, it was assumed that any dynamic loading applied by the IHS 
would be transferred through the center portion of the plate to the preload 
screw without affecting the stresses in the plate produced by the HSHX loading 
around the outer edge of the plate. 


76 






















































ffesIgn^VoIablT^s^V “IhVma^T* I" tab ' 6 ?8 ’ *' S ° sh °“" »• the 

between the • 1^1 .?*! ined as the difference 

~r„. «%.“£ ’ZZZFfzz S 

»» is.tr .r«" ■ srsyisr air ■» - 
;;;:rr£ ■=£ HH-3": 1 sat 

response at the HSHX mount points with ih e V ^ S "k S < T°" ducted + ° determine the 
inputs from MIL-STD-810B applied to the HSA mounts ^ b °° S+ sinusoidal 

were 60 Hz and 118 Hz for the HSA and ?61 frequencies of interest 

sinusoidal input at these friueniii! u T ' 2 f ° r +he HSHX ‘ The 
The load path from the HSA ^unts lo th^ HS^nl^ "V? 70 ? 108 95 ± 3 * 5 9. 
joints that can provide damping. The specif fe^dal^n I* i nvolves several 
thought to be conservative and was used for the ? ? factor of 5 percent was 

HSHX mount points was determined for both axial and tat S ’ i™* response at +he 
results ot this analysis are sunmarlted In tible 29. « c,t atlon. The 


TABLE 29 

HSHX MOUNT RESPONSE TO MIL-bTD-810B INPUTS 


Excitation 
d i rect ion 

Frequency, 

Hz, 

Input to 
HSA, g 

Response at 
HSHX mounts, 3 

9 

Amp 1 i f i cation 
factor 

Axial 

(Z) 

Latera 1 
(X) 

60 

118 

261 

342 

60 

1 18 
251 
342 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

5.6 
15. 1 
1.8 
1.3 

3.6 
4.0 

6.7 

7.8 

1.6 

4.3 
0.5 
0.4 

1.03 
1.15 
1.92 
2.22 

— L 



a With 5 percent damping 


The maximum amplification between +h« uc* 
points is 4.3 at the 118-Hz axial resonlnt T m ° UntS and the HSHX attachment 
amplification at the HSHX resonant frequenc i e^ i T???' HSA * The maximum 

factors result in a 15.1-g maximum loS at the HSA-HSHX ill* 5 ? ^ " ! lc3t ion 
•s well below the 50-q equivalent docil . 2 "SA-HSHX interface. This load 

HSHX mount brackets and the HSA support' p,"e (see fTtuU?*'**'* ° f ^ 
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Type of 
stress 



Applied stress 

Margin 
of safety 

Component 

Mater i a 1 

Design criterion 

MPa 

psl 

Pin , 4.78 mm 
(0.188 in.) 

Bending 

Waspaloy 

Endurance limit 
310.3 MPa (45 000 psl) 

279.6 

49 560 

0.11 

Bracket, 4.06 mm 
(0.16 in.) 

Bending 

Hastel loy X 

Bearing yield 

335.4 MPa (48 640 psi) 

139.1 

20 180 

1.41 

Bracket weld 
2.79 mm (0. 11 in. ) 

Bend i ng 

Hastel loy X 

Endurance 1 imit 3 
135. 1 MPa ( 19 600 psi ) 

67.1 

9 730 

1.01 

Braze at base 

Ten s i on 

Palniro 7 

Endurance 1 imlt a 
34.5 MPa (5 000 psi) 

P.9 

1 730 

1.89 

Aft support plate 

Bend i ng 

Poco graphite 
AXF-5Q 

Endurance limit 3 
77.6 MPa (11 250 psi) 

56.1 

8 140 

0.38 

Aft support plate 

Bearing 

Poco graphite 
AXF-5Q 

Bearing yield* 5 
103.4 MPa (15 000 psi) 

47.7 

6 920 

1.17 

Aft support plate 

Shear 

Poco graphite 
AFX-5Q 

Shear yield 0 

51.7 MPa (7 500 psl) 

4.8 

690 

High 


a Estimated endurance limit = 0.75 x yield stress. 
Estimated bearing yield = yield stress. 


c Es+imated shear yield = 0.5 x yield stress 
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system design 


Design Description 

The design effort in f“atiws,^ «^ ny 

Tb t e S si m llar r “ a 'i?s e p'-«ent design as specified by the Genera 

The Key recommended HSA 'incase" the natural tr ° q “^ cy °* 

r e *vit th tbe 9 «Hr c :°un u t r ;:,i^ sir. : rssJt™ 

sKet* of the end enclosure Is Shown f present I y^composed of 60 layers 

navi been made tor mod 1 1 ica t.on * « •* V.ever . that changes will -e ° 

°* rlTord^T^t-s^itred P?CS temperature limits for E CS opera 

The HSHX is a single-tube [tosteHoy * h “' InW iaM y h'ckness. 

s 

ducts* ar^short^ Haste, « ZSZZZv* 

r as i°;I 0 o°tu n r;l rt'laThSnJ'th. coil near the elbows to prevent excess, ve 
bending in the tubing. ce between turns of 

The HSHX coiled such that there ;s a ncm (0 .„ in., wide 

. ,n rv?R in ) A continuous Haste i loy » K . h rni i The brazed strip 

n 71 mm (0.028 in.j. ^ . hraze (j to tne con. 11,0 “ 

" coi 1 • A 5ke,ch of the braz c 

assembly is shown in tig. 31. 

♦in ^2 The forward mount 

ing ^;^"°{^hets-rat-n r :::r 

ec.W- - ° 0E and are 

other minor changes 
discussed in ref. 10* 

Photographs o, the HSHX assemb, y. Including the mounts, are shown ,n 

,i9 ‘ tlrB hSA with emphasis on the HSHX, has been 

A design layout tor the e t HS ^ of th |s report, 

prepared and is included as t,g. 5“ 
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F i quire 31. — Brazed Coil Ass 




Figure 32.— Aft HSHX Mounts 


PAG! IS 

OF POO* QUALITY 






PfcC »<*.2 


.-T- 


P60S9^-1 


i^nre 33. HSHX Assembly. 


DIOGINAS b 

M eoor QUAUTY 


,J *Ai 


W nr, Pa 



v 






1 


l 




HSHX Manufacturing Procedure 


The welded and drawn Hastel loy X tubing should be Inspected visually and 
radiographically to ensure that there are no significant material flaws. Trace- 
ability and lot control should be established. The tube Is then formed Into a 
coll with particular attention paid to avoiding flat spots and unacceptable 
wall thinning due to bending. The transition elbows are then welded to the 
coll. Visual and radiographic Inspections should be conducted to control and 
verify weld integrity. A helium leak test should be performed at ambient con- 
ditions to ensure integrity of the entire tube. 

Assembly of the HSHX Is expected to utilize two braze cycles. For the 
first braze, the continuous coll strip and the reinforcement strips are 
brazed to the tube using the gold-based brazing alloy, Palniro I. During braz- 
ing of the coll, special fixtures should be employed to hold proper coll dimen- 
sions and maintain the continuous strip In its proper location. The brazing 
is performed In a vacuum furnace at around 1132*C (2070*F). The mounting 
brackets are joined to the coil in the second braze cycle, using another gold- 
based braze alloy, Palniro 7, which has a lower melting point than Palniro 1. 
Prior to this second braze, the bracket assembly Is contoured and positioned 
against the coil and the brackets tack welded to ensure proper location during 
the braze cycle. The braze Is performed at around 1052*C (1425*F). 

Following the brazing operations, detailed radiographic inspections should 
be performed to verify braze quality. After inspection, installation of the 
mounting pins and final machining operations are completed. The pins are 
installed after brazing to avoid any changes in the Waspaloy heat treatment. 

The mounting pins are tack welded to the brackets. 

After all manufacturing processes are complete, acceptance tests should be 
performed to verify unit integrity and operation. These should consist of an 
isothermal pressure drop test and a proof pressure test. An additional helium 
leak test should be performed at ambient conditions. 

HSHX assembl ies of this design have been fabricated for use in the BIPS 
ground demonstration system (GDS). The basic manufacturing feasibility has, 
therefore, been veriflod. The coi l-formli,g process required development 
because of the combination of tight coil radius and multiple-turn configura- 
tion. Satisfactory parts were produced by repeatable procedures and are used 
in The GDS HSHX assembl Ies. 

Heat Source Assembly Weight 


The expected HSA weights are listed in table 30. The weights are for 
the reference design HSHX. The recommended modifications to tta HSA end 
enclosures are included. Excluding the IHS, the heaviest components, and 
hence those with the greatest weight savings potential, are the HSHX and the 
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TABLE 30 

hsa weight SUMMARY 


HSHX 

Housing assembly 
Insu I at ion cy I inder 
Insulation ends 
End enclosures 
Preload screws 
Domes 

Miscellaneous hardware 
Subtotal 
IHS 

Total 


—ha.. 

Jb 

10.07 

22.20 

3.23 

7.11 

7.03 

15.49 

1.41 

3. 10 

2.26 

4.98 

0.29 

0.65 

1.21 

2.66 

3. 18 

7. Qi 

28.68 

63.20 

18.82 

4 1.5Q 

47.50 

104.70 


~ saw* 

J ° '■ 5 lb > for each HSA. 
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R&OWENDED WORK 


Future work in the following areas is considered to be of importance to 
the final design of the HSA and the verification of system performance. 

Analytical Work 

Investigation of theBIPS startup transient, HSA fueling, and ECS 
operation for a HSHX with a surface emissivlty of 0,25 (untreated) . — Th i s 
emissivity has been shown to be advantageous for maintenance of the PICS 
lower temperature requirement for normal operation and ACS deployment. 

Tradeoff study to optimize the HSIS for normal operat ion .--HS IS weight 
(number of fo< I layers) should oe varied to determine the effect on HSA 
heat loss. Tradeoff factors would have to be supplied by the 3IPS program 
office. The effect of varying the number of foil layers on the other nodes 
of system operation must be considered. 

Investigation of HSA cooling prior to RIPS startup .— The latest conception 
of the 3 1 PS f I ight system el iminates the cocoon surrounding most of the HSA 
and BIPS components. This results in a surface cooling requirement for the 
entire HSA during launch pad ACS operation to meet the exposed surface 
temperature limitation. In addition, a maximum heat dump limitation to the 
spacecraft may require cooling of the HSA both on the launch pad and in space 
prior to startup. The various cool ing requirements need to be Integrated and 
a heat exchanger designed to provide this cooling (probably a coil wrapped 
around the outside of the HSA case). 

ECS operat I or . — The tnermal analysis conducted during this study 
indicates that the PICS temperature limitation for the post-meltdown steady- 
state condition cannot be met by the present ECS. Lower melting point toil 
materials should be investigated in an effort to melt more of the foils. 
Combinations of higher and lower melting point foils might be used and the 
overall number of layers could be reduced. A detailed effort should include 
the effects of foil sublimation, which may reduce the number of foil layers 
remaining after an extended period of time. Revis'on of the thermal model 
to .nclude a single mode for each foil layer may be worth the resulting cost 
in computer time due to the increase in accuracy. 

If the PICS tenperature goal cannot be met using the foil meltdown ECS, 
alternate ECS concept s m ight be investigated. This might be necessary, since 
essentially all of the foils must be reroved to reach the steady-state PICS 
temperature limit. Most of the alternate systems involve the bulk removal of 
the foil insulation. Those methods have the advantage of decoupling the HSIS 
and the ECS, thus allowing greater flexibility in HSA optimization. Foil 
removal by meltdown would act as a backup for any alternate ECS. 

Investigation of other operating conditions . — A matrix of potential 
operating conditions for the BIPS flight system is available. Some of these 
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operating conditions may warrant complete thermal analysis for prediction of 
PICS temperatures or as inputs to a structural analysis. 



Increase in turbine inlet temperature by shielding of fhe HSHX .— By 
shielding or insulating the hot end of the HSHX, more of the heat released 
by the fuel is absorbed in the colder regions of the heat exchanger. The 
total heat absorbed by the working fluid will be essentially unchanged, 
since the HSIS effectively determines the HSA heat losses. Thus, the HSHX 
outlet temperature can be increased with no increase in maximum tube wall 
temperature. The shielding acts as a re-radiating surface for the heat 
emitted by the fuel; the re-radiated energy is eventually absorbed in the 
colder HSHX region. 



Various shielding schemes should be investigated for both physical 
arrangement and fluid temperature increase. The optimum situation tor fluid 
temperature would be to distribute the heat so the entire heat exchanger 
surface facing the I HS is at a uniform temperature. The heat flux at the 
hot aid of the HSHX would be very small, so that the fluid outlet temperature 
could be very close to the maximum wall temperature. This uniform wall 
temperature condition might be approached by telescoping the insulation 
around the HSHX, with few layers of insulation in the inlet (cold) region 
and a greater number of layers in the outlet region. 

The shielding may have some unfavorable effects on the system, such as 
increasing the PICS temperature and increasing the circumferential temperature 
gradient in the cold end of the HSHX. These potential effects would have to 
be carefully investigated. 

Stru c tural ana I ysi s . --Further analytical work would appear ro be necessary 
if a more accurate assessment of LCF life of the HSHX is desired. More 
detailed characterizations of the localized stresses in the HSHX coil and mount 
interfaces would be desirable. The thermal stresses produced during the HSA 
fueling transient have not been evaluated for their effect on fatigue life. 

In addition, structural analyses could be performed for other BIPS flight 

system operating modes. -I 

The fatigue life of the remainder of the HSA has not been investigated 
in detail. An area of particular interest Is the response of the graphite . I 

support plates during the HSA fueling. Imperfect thermal contact between 
the plaTes and the I HS may cause significant thermal gradients and stresses 
in the support plate. An investigation of the proposed flight system 
operating matrix may indicate other modes that should be studied for effect 
on HSA LCF life. 

Experimental Effort 

Stability of PICS and HSHX emisslvities .— The iridium PICS is surface- 
treated to increase its emissivity. The stability of the surface over long 
periods of time should be investigated since the PICS temperature is very 
strongly dependent upon its emissivity. 
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Similarly, the HSHX temperatures are dependent upon the heat exchanger 
surface emissivity. The reference design HSHX calls for an emissivlty of 0.8; 
considerable surface treatment may be necessary to attain this emissivity. 

To ensure that the HSHX temperatures do not exceed the requirements during 
the BIPS operating life, the stability of the surface must be ascertained. 


HSIS heat transfer correlations . — The representation of heat transfer 
through the insulation could be improved. Data are especially needed for 
heat transfer through foils where the cold end is well above ambient tempera- 
ture. An improved correlation also is necessary to account for the edge losses 
through the insulation joints. 


PICS temperature measurement . — Direct measurement of the PICS temperature 
Is necessary to confirm the thermal model of the IHS. Since the model Is the 
same for steady-state and transient (including ECS) conditions, steady-state 
experimental results would suffice to confirm the PICS temperature prediction 
for all operating conditions. The IHS could be used by itself, with carefully 
defined boundary conditions and axial measurements of the PICS and ablation 
sleeve temperatures. 


Structural verification . — In addition to the normal acceptance- type tests, 
specific testing should be performed to confirm the dynamic response and 
capability of the HSA and HSHX. In particular, the amplification of dynamic 
loads at the HSHX mount points is an important area of investigation. 


1)1 
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SUfflARY OF RESULTS 


The following Itemized remarks summarize the results of the study: 

(1) A Hastelloy X single- tube helical coll was selected as the design 
tor the heat source heat exchanger (HSHX) for the Brayton Isotope 
power system (BIPS). The material and configuration emphasizes 
reliability and tabr Icab 1 1 1 ty. 

(2) The feasibility of fabricating this design has been demonstrated. 

A heat exchanger of this configuration has been manufactured for 
use In the BIPS ground demonstration system under a TOE contract. 

(3) The single-tube hel icai-col I design has been optimized with respect 
to weight, pressure drop, and working fluid outlet temperature. 

The optimized heat exchanger is a 10-turn col I made of a 3.%-cm 
(1.56- in.) 00 and 3.703-on (1.458- in.) 10 tube. The surface should 
be treated to an emissivlty of 0.8. 

(4) The HSHX meets the steady-state performance requ i rements. The 
temperature requirements for the Iridium post- impact containment 
shell (PICS) are not quite met. The FICS temperatures are within 
the required range for an HSHX surface emlsslvity of 0.25, the 
nominal value for untreated Hastelloy X. 

(5) The auxil iary cool ing system requirements can be met using neon 
as the fill gas. The neon fills the radlatio.i gaps In tho hoat 
source assembly (HSA) and allows the heat to be transferred out of 
the HSA by conduction. Additional cooling of tho HSA Is required 
to meet a launch pad exposed surface temperature limitation. The 
PICS temperature requirements are not quite met for the reference 
HSHX, but are within the required range for an HSHX of 0.25 surface 
emlssl vlty. 

(6) The latest expected BIPS startup transient was used to predict 
temperatures, thermal stresses, and low-cycle fatigue life. The 
PICS temperatures are acceptable during the startup. Tho predicted 
fatigue life of the HSHX is 570 cycles, well over tho system 
requirement of 100 cycles. 

(7) Tho emergency cool Ing system performance Is not sufficient to meet 
the PICS temperature requirements. Although tho maximum allowable 
PICS temperature Is not exceeded during tho foil moltdown, an 
insufficient number of foils Is melted to lower tho ultimate post- 
meltdown PICS tomperaturo to Its required value. Analysis indicates 
that essentially all of tho Insulation foils must be removed to 
attain the required post-mel tc.own PICS temperature. 

(8) With some minor modifications to tlio oxlstlng HSA structure, natural 
frequency and dynamic loading requirements can be met. 
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CONCLUDING REMARKS I 


This final report has presented the results of a study conducted by the 
AIResearch Manufacturing Company of California on the preliminary design of 
an alternate heat source assembly (HSA) for a Brayton Isotope power system 
(BIPS). In particular, the HSA is Intended for use on the BIPS now under 
development by the AIResearch Manufacturing Company of Arizona under a 
Department of Energy contract. 

The heat source heat exchanger (HSHX) design for the HSA has been selected 
and optimized. The auxiliary cooling system has been investigated, and neon 
recommended for the HSA fill gas. Low-cycle fatigue life requirements are 
easily attained for the HSHX. Natural frequency and dynamic loading require- 
ments can be met for the HSA. 

The major problem area remaining for the HSA Is the emergency cooling 
system (ECS) performance. An Insufficient number of insulation foils is melted 
to lower the temperature of the iridium post-impact containment shell to an 
acceptable value. 

Reconmendations for future analytical and experimental work, with 
emphasis on the ECS problem, have been made. A layout drawing has been pre- 
pared for the HSA. 
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Program Capabilities 

, s a j Research thermal analysis program H09I0. has the 
This program, known as 

following capabil Transient heat transfer calcula- 

^SSsSSSSssss rw~ 

^-a r «e.t. heat trans jer :c alculatl^ : ;-St«^ 
c ulat^ons^are"*based "on a teihnlgue Involve, "acce erated 

ef a! t tut ion «lth monotonlc deceleration - n !" s ore s|0 . t0 

also provides tor at.y ° 

tion.^or radiation heat transfer. , rans(er ^ ^ 

5«^31»sSKr.ts 

the thermal conduct iv y function of temperature, 

tied as a constant value or as a tunc 

Convection ±eat^ ^ Th'^solTT' 

,0 th^^hy sp r l.yih9 r » s Sond t elemen arM ^ the 

I a %\:\n«.':?th 0 cony^lon from the solid surface to the fluid.. 

0 the conuect,on r , transfer 

.1 r,f the "film" temperature In a . nc i en t tor both open and 

tion of The < 4i _ he«t transfer coetticien nrocjram 

^jrr=s--s=!Jf as*---. — - 
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. tree or «c— rotating 

•‘windage” heat generation. In met od 9 from rQW of ho | es onto a 

ment heat transfer coeffic ®" * t ° pr ^gr«n calculates convection heat transfer 
concave surface. In method nine P ^ transfer "fin effective- 

coefficients for flow In a d ^t including the^ utn|zes tables of 

ness" of extended surfaces "> £ f Reynolds nunber to the program. These 

Colburn J-factors Input r L,d ducts, square ducts, rectangu- 

tables may be generated for flu es dimpled tubes, and curved ducts, 

lar ducts, triagular ducts, annu P ! Jn ^ banks# plate-fin surfaces, 

They also may be generated for flu crossed rod matrix surfaces, and cor- 

pin-fln surfaces, screen *®rt effects on heat transfer may be applied 

rugated ceramic surface 5 . En tra location. Four techniques tor 

using the appropr »*• mu. ^i^^Id^ndent fluid properties are 

available 0 In^ tna* program" The appropriate fluid properties nay be Input In 
tabular form as a function of temperature. 

Cd.atlnn neat transfer caJ al^.-^dWI^^t ^anster^lsjnput^o 

the program by specifying a solid * **•" £ croS s-sectional area for radiation 
surrounding element number b Y rad ' a J fl , th from the center of the solid 
from the solid element, and ^ ond “ ^Jalso includes the important 

element to the radiation surface. * surface with radiation from the sur- 

combination of conduction to ^e f ed latFonsur face ated by metho ds 

irz Misii a computer pro9rw 

such as CONF AC I I. 

_LELULL®J ^©rn£or_a+_ur©^ 

and thermal co nduct . v i ty spec | cat »o • _ flnd thermal conductivity may be 

conditions, heat input, thermal capa ce, Qf elemen ts that are iden- 

specified for each mdividu o.i rulation the initial temperature, the 

tical • in transient heat transfer ca ' ^'at ion, the^l f ^ ^ ^ +herrna | ^ 

heat Input, the dwsity, the F £ r 9l(men , s with negligible thermal 

ductivity of each element is specif t mg be left b | a nk to 

capacitance the density, volume, and specific d ^/calculations, the Initial 
increase the calculation time P* . conductivity of each element are 

temperature, the heat input • , 2™ aes£c i f I ed by having a negative 

specified. The boundary cond It ion el enents a maintained at a 

value for the density times the vol umo. verS us time function 

constant temperature or may be s P®c' f >® d * J P be specified as a boundary 
from an input table. Any element n the network may^ J ^ be con . 

condition (constant temperature) element 3 Jadi^tion. The heat input tor each 
nected to it by conduction, convect ion or red let .on. ^ Qf + ime ln a 

element may be input as *® ro ' t(OTDera ture or another specified element tern- 

table, as a function of 15 ° ••wlndaqe" heat generation calculations, 

perature. specified from the ?" 8 ' “hJ* 9 heat generation calculation com- 

or calculated from the ball 0 specific heat and thermal conducti- 

v rtror^r.rrnt C : n ay be be S s P p P ei:rl;d SVSifnt or as a function Of 

temperature In tables. 
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Fluid stream heat transfer and pressure drop calculations . — Fluij strean 
elements may be input with heat transfer to them by conduction, convection, or 
radiation. Fluid stream heat transfer calculations have provisions for pre- 
venting the outlet fluid temperature from "overshooting" the surrounding sur- 
face temperatures, a thermodynamic impossibility. The steady-state fluid 
stream calculations are based on thermal capacity rate calculations, while 
transient fluid stream calculations may be based on the thermal capacitance of 
each element moving in the fluid stream to simulate "lag" conditions. The 
energy Input of rotational flow may also be added to the fluid stream. 


Both steady-state compressible and incompressible fluid stream pressure 
drops may be calculated by the program. The pressure drop calculations include 
the effects of heat addition, area change, fluid friction, rotational flow, and 
flow addition or removal. Total head losses due to valves, bends, sharp con- 
tractions or expansions, and orifices may be Included at the inlet and exit to 
each fluid stream. 


A complete fluid stream network may be simulated with streams branching 
from previous streams and mixing to form new streams or even returning to a 
previous stream in the network. The fluid flow rate may be input as a con- 
stant, as a function of time, from a table, or as a function of specified 
element temperature. 


Program Output 

The following outputs may be obtained: 

(1) Each element temperature, heat input, and thermal conductivity for 
steady-state calculations is printed out. Each element temperature, 
heat input, weight, specific heat, and thermal conductivity for 
each specified printing time period in transient heat transfer 
calculations is printed out. 

(2) The fluid stream inlet temperature and the outlet temperature, the 
fluid stream flow rate, the fluid density, and the internal fluid 
heat generation for each section of each fluid stream are printed out. 

(3) The "free stream" temperature, the film discharge temperature, and 
the effective film temperature at each location specified are printed 
out . 

(4) The printing of the thermal resistance values for conduction, the 
thermal resistance values and heat transfer coefficients for convec- 
tion, and the thermal resistance values and effective heat transfer 
coefficients for radiation may be included or deleted as specified. 

(b) The fluid stream pressure drop calculations and printout may be 

deleted if specified. When Included, the total and static pressure, 
the Reynolds number, the friction factor, and the Mach number for 
compressible flow are printed for each element In each fluid stream. 
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Typical Applications 


Typical applications include: 

(1) Both passive and active electronic cooling system analysis and design 
with or without heaters or cooling flow controllers. 

(2) Thermal analysis and design calculations tor ambient-cooled, forced- 
air-cooled, gas-cooled, or liquid-cooled ac or dc motors, generators, 
and alternators. 

(3) Thermal analysis and design calculations tor pumps, fans, and 
compressors, including the bearing temperatures, and the analysis 
of the motors or turbines driving them. 

(4) Thermal analysis and design calculations of gas turbine engines 
including the axial-flow and radial-flow compressors and turbines, 
the bearings and seals, the anti-icing system, the lubrication cool- 
ing system, the fuel sypply system, and the accessory area cooling 
system. Also, the thermal analysis of cooled and uncooled turbine 
bl ades. 

(3) Transient and steady-state thermal analyses of heat exchangers 
including air-oil coolers, fuel-oil coolers, recuperators, 
rotary regenerators, cryogenic heat exchangers, pool boiling 
heat exchangers, condensers, periodic flow regenerators, and 
heat exchangers with more than two fluid streams. The calcu- 
lations may Include the effects of axial conduction, fluid 
bypassing, perfectly mixed or unmlxed fluids, variation of 
fluid properties through the heat exchanger, condensation of 
moisture from the air or "wet" heat transfer, and the effect 
of the variation of f I uid-to-wal 1 temperature difference on 
local heat transfer coefficients tor boiling and condensing. 
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U-FLOW MANIFOLD DESIGN COMPUTER PROGRAM, XIOOO 


-err sxrrsr is 4r r-Kr,?"'.:- 

s-JSiSfs sna :;-,k -• 

calcu!ater eS The\Jh^ ribU !!°?^! 0n9 th ® C ° re ,en9th of the 9 iven manifold Is 

«“ ir-rKi^s^-si ss^v r ,le - 

motion. Either one- or two-dimensional solutions can be specified Rnth 

in"?l?fo,d n s COmPr : SSib, r Ca,CUl0t, - S «" be^erformedh^Fr let iona^^osses 

in the man folds can be included by the input of appropriate friction fJw 
versus Reynolds number tables for each manifold. 

Program input comprises the following: 

Fluid properties 


Manifold friction factor properties 

Specification of options and convergence criteria 

Definition of one manifold flow area and hydraulic radius 
as a function of manifold i ength 

Mass flow rate, inlet and outlet temperatures and pressures 
Program outputs comprise: 


Design manifold definition 
as a function of manifold 


- flow area and hydraulic radius 
I ength 


Difference in manifold pressure profiles as a 
manifold length 


function of 


Maximum heat exchanger core pressure drop variation for the 
design solution 


Manifold | ressure losses 
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ANSYS (ENGINEERING ANALYSIS SYSTEM) COMPUTER PROGRAM 
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